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I. Sununary of Work 


The original objective of the work was to obtain radiometric measure- 
ments of cloud Influence onUV-Band define mathematical models of the In- 
fluence so as to lay the groundwork for later construction of the global 
UV-B climatology from satellite-determined ozone data. We believe we have 
made considerable progress towards these objectives as described In the 
attached papers. 

The related radiometric work first undertaken by the group used all-sky 
photographs In conjunction with an Eppley radiometer and a Robertson-Berger 
meter. This particular study. Initiated under the auspices of the CIAP 
program of the Department of Transportation but completed during the period 
of the grant, Is described In Attachment 1. This work represents our best 
.attempts at the broken cloud problem. More refined measurements comparing 
UV-B radiation with total solar radiation were carried out somewhat later 
in conjunction with the theoretical analyses given in Attachment 2. In 
Attachment 3 the cloudy case is referred to the cloudless sky irradiance 
and convenient "transmission" ratios are given. Attachment 4 summarizes 
an approach to the inversion of scattering data which should eventually 
facilitate the final objective of this work. Attachment 5 presents an 
improved characterization of the UV-B radiation from a cloudless sky. Such 
a result is not only of use in itself but also should facilitate further 
cloudy sky analyses. Attachment 6 extends rbise results to include actinic 
flux. 

These studies along with studies undertaken by the writer in connection 
with his work as a member of the Nimbus-7 SBUV/TOMS experimental team 
have laid the groundwork for the later construction of a global UV-B 
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cllmatology from satellite-determined ozone data. We hope eventually to 
obtain funds to attempt to satisfy this final objective with the aid of 
Nlmbus-7 UV data when It becomes available. 

Not the least contribution of our program of UV studies has been our 
output of trained manpower, many of whom received their first exposure to 
atmospheric radiation research in our group. Dr. An Ti Chai is now at the 
NASA Lewis Research Center wottcing on solar radiation and solar cells. 

Dr. Tsan Mo Is now working at Goddard Space Flight Center. Dr. Janusz 
Borkowskl has returned to the Institute of Geophysics at the University of 
Warsaw, and continues to be Involved In stratospheric problems. Dr. James 
Splnhlme and Dr. Richard McPeters who received his Ph.D. degree during this 
period are also at NASA Goddard Space Flight Center working on atmospheric 
radiation problems. Dr. Kenneth Klenk now with System and Applied Sciences 
heads a project on the processing of Nimbus-4 and Nimbus-7 SBUV/TOMS data. 

Mr. Kenneth Cross is now working in industry as a program analyst. Mr. L. R. 
Smith, an undergraduate completing his B.S. degree, has undertaken a senior 
project in atmospheric sciences. Mr. David D. Doda and Dr. Dale Brabham 
shofld be completing their programs with us by next summer. Mr. Paul 
Schippnick has recently been initiated into theoretical atmospheric radiation 
studies. We hope to add an experimental graduate student with industrial 
experience in the next two months. 

This Final Report of our study on the "Influence of Clouds on UV-B 
Penetration to the Earth *s Surface" for NASA Contract NAS9-15114 covers the 
entire contract period. This was originally intended to be from 
September 1, 1976, th?.ough August 31, 1978. However, because of scientific 
manpower shortages the program was stretched out and extended to September 30, 
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1979) with no additional funds beyond the original $50,000 award. The 
stretch-out was also facilitated by the fact that during the period of 
this contract, the principal Investigator was also supported by other 
grants and contracts, some of which overlapped the coverage of this contract 
In various ways. Section II contains a list of publications of the 
principal investigator since 1976. The Items encircled represent works 
closely related to the topic of this grant. The Items encircled and 
surrounded by hexagons are directly related to this contract. Copies of 
these publications are attached. Also attached Is the Summary Chapter of 
a forthcoming book on a national problem of potential interest to NASA. 

The book concerns a topic which potentially has greater climatic and health 
consequences than the UV-ozone problem which has been the focus of this 
study. 

We trust that this Final Report fulfills the requirements of this 
contract and look forward to future collaborations with NASA-Houston. 
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Ratio in the Middle Ultraviolet" with Kenneth F. Klenk, Proceed- 
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^42^ "Spectral Sunphotometer Using a Compact Spectrometer", 
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"Florida's Air Quality, Present and’ Future," with D. F.. Rio, and 
R. A. Hcdlnger, Florida Scientist, 3, 182-190, (1978). 

"Total Ozone Determination by Spectroradiometry in the Middle UV," 
with L. M. Garrison and 0. U. Doda, Appl. Optics, 18 , 6, 850-855, 
(1979). 

"Remote Sensing of Atmospheric Aerosols," University of Florida 
Press, with D. A. Lundgren, Ed., 635-650, (1977). 
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Control," with E. T. Loehman, S. V. Berg, A. A. Arroyo, R. A. Hedinger, 

J. M. Schwartz, M. E. Shaw, R. W. Fahien, V. H. De, R. P, Fishe, 

D. E. Rio, and W. F. Rossley, J. Environmental Economics and Manage- 
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"Electron Impact Excitation of the Beryllium Isoelectronlc Sequence," 
with P. S. Ganas, Physical Review, 19 , 6, 2197-2205, (1979). 

"Remote Sensing of Ozone in the Middle Ultraviolet," with J. D. 

Talman, to be published in Proceedings of Workshop on Interpretation 
of Remotely Sensed Data, Williamsburg, Virginia,' May 23-25, 1979, 
Academic Press, (1980),, 

"Electron Impact on Atmospheric Gases 3. Spatial Yield Spectra for 
N_," with C. H. Jackman, J. of Geophysical Research, 84 , A6, 

2715-2724, (1979). 

"Comparative Studies of Atomic Independent Particle Potentials," 
with J. D. Talman and P. S. Ganas, International Journal of Quantum 
Chemistry, in press. 

"Pollution Dispersal Modeling and' Regional Public Policy Impli- 
cations," with E. T. Loehman, S. V. Berg, R. W. Fahien, M. E. Shaw, 

M. J. Jaeger, H. Wittig, A. A. Arroyo, J. M. Schwartz, R. A. 

Hedinger, V- H. De, D. E. Rio, T. J. Buckley, R. P. Fishe, W. F. 
Rossley, and D. Trimble, to be published in Proceedings of Sym- 
posium held in Williamsburg, April 17-19, 1979, on Environmental 
and Climatic Impact of Coal Utilization, Academic Press, (1980). 

"Microplume Model of Spatial Yield Spectra," with R. P. Singhal, 
Geophysical Research Letters, 7, 625-628, (1979). 

"Factor of Safety Method Application to Air and Noise Pollution, 
with T. J. Buckley, D. E. Rio, A. MacEnchern and R. Makarewicz, 
to be published in Atmospheric Environment. 

"Improved Analytic Characterization of Ultraviolet Skylight," with 

K. R. Cross and L. A. Smith, to be. published in Photochemistry and 
Photobiology. 

"Correlation Between Noise and Air Pollution," with R. Makarewicz, 
to be published in Sound and Vibration. 



"Vertical Transport of N0_ Through the Tropopause," with J. 

Borkowskl, published in ACTA GEOPHYSICA POIGNICA, 26, 4, (1978). 

"Aerosol Effects on Atmospheric Radiation," with P. F, Schlppnick, 
to be published in Proceedings of Workshop on Atmospheric Aerosols: 
Their Formation, Optical Properties and Effects, held at Baltimore, 
Maryland, November 6-8, 1979, Academic Press, (1980) . 

"Spatial Aspects of Low and Medium Energy Electron Degrddatlon in 
N 2 »" with R. P. Slnghal and *C. H. Jackman, to be published In Jour. 
Geophysical Research, (1980). 

Coal Burning Issues , edited by A, E. S. Green, to be published by 
University Presses of Florida, 400 pp (1980). 

"Atmosphere Modification," with K, E* Taylor and W. L. Chameldes, 
Chapter 11 of Coal Burning Issues (see Number 268) , 203-230 (1980) . 

"Quantitative Public Policy Assessments," with D, E. Rio, Chapter 15 ^ 
of Coal Burning Issues (see Number 268) , 303-330 (1980) . 



III. Recent Contracts and Grants 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

1. Aurora and Alrglow Phenomena in Planetary Atmospheres, NGL-10-005-008. 

Fundamental Investigations of alrglow and auroral phenomena in 
planetary atmospheres. $50K, 5/1/79-4/30/80. 

2. Science support for the Nimbus-7 Sensor, NAS-5-22908. 

As member of the Nimbus-7 Experiment Team (NET) for the Solar and 
Backscattered Ultraviolet and Total Ozone Mapping System (SBUV/TOMS) 
the P.l. assists iia the determination of solar ultraviolet irradiance, 
ozone profiles, and global maps of total ozone. We provide experi- 
mental measurements of the ultraviolet levels reaching the ground, 

UV reflectivity measurements and total ozone values for correlation 
with the satellite measurements. $208K, 12/19/75-10/1/80. 

3. Influences of Clouds on UV-B Penetration to the Earth's Surface, NAS-9-15114. 

Define mathematical model of UV-B irradiance reaching the ground and 
the influence of clouds. $50K, 9/1/76-9/30/79 - expired. 

NATIONAL SCIENCE FOUNDATION 

1. Atmospheric Radiative Transfer and Atomic Processes, NSF-ATM-75-21962. 

Theoretical and experimental atmospheric radiative transfer studies 
and atomic calculations for aeronomical applications. $39K, 5/26/78- 
10/31/79 - expired. 

U. S. DEPARTMENT OF ENERGY GRANT 


1. Biophysical Studies Related to Energy Generation, DE-AS05-76EV03798. 

Investigation of the physical stage of the interaction of the radiation 
with the inert matter and living organisms and on dose-response re- 
lationships especially related to energy generation. $50. 3K, 6/1/79-5/31/80. 

ENVIRONMENTAL PROTECTION AGENCY 

1. The Impact of Stratospheric Ozone Depletion upon Tropospheric Ultraviolet, 
Photochemistry and Smog, EPA-R806373010. (With W. L. Chameides) $50K, 
5/1/79-4/30/80. 

Radiative transfer calculations and modeling studies are used to determine 
the impact of stratospheric ozone depletions upon the tropospheric UV 
radiation field and photolysis rates, tropospheric photochemistry and smog 
formation. 

UNIVERSITY OF FLORIDA (GATORADE) . ICAAS 

1. Comprehensive Coal Burning Issues Study (with 30 U.F. faculty members). 

ICAAS is conducting a detailed examination of energy and environmental 
issues facing the state and the country during the next quarter century, 
particularaly as it relates to coal burning and its atmospheric impacts. 
$150K, 7/1/79-9/30/80. 

FLORIDA BOARD OF REGENTS 

1. Impact of Increased Coal Use in Florida, ICAAS, STAR grant 79-064. 

$67K, 7/1/79-9/30/80 (with 15 U.F. faculty members). 

ICAAS is studying the major economic, social, environmental and techno- 
logical issues involved in the anticipated increased use of coal in Florida. 
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CALCULATION OF THE RELATIVE INFLUENCE OF 
CLOUD LAYERS ON RECEIVED ULTRAVIOLET 
AND INTEGRATED SOLAR RADIATION* 

J. D. Spinhirne and A. E. S. Green 
U niversity of Florida, Gainesville, Florida 32611, U.S.A. 

(First received 1 March 1978 and in final form 6 June 1978) 

I Abstract - Discrete ordinate radiative transfer calculations are used to examine the relative influence of cloud 

layers on the received ultraviolet flux and the received solar energy for a plane parallel atmosphere. The 
wavelength dependence for the cloud influence on the integrated solar radiation is primarily a function of the 
interaction of the cloud layer with the surrounding atmosphere and the underlying surface. It is found that for 
u.v. wavelengths greater than 300 nm, the ratio of u.v. to solar energy flux transmission through the 
atmosphere is insensitive to changes of cloud height, cloud scattering parameters and surface albedo (for 
realistic values), but is dependent on cloud thickness. For wavelengths approaching and less than 300 nm, 
absorption by tropospheric ozone results in sensitivity to cloud height and surface albedo. 


1. INTRODUCTION 

Ultraviolet (u.v.) radiation is known to influence many 
biological processes and it is also known that the 
earth’s ozone layer controls the amount of u.v. ra- 
diation which reaches the ground. In response to the 
great concern that anthropogenic modifications of the 
ozone layer might have significant health and agricul- 
tural impacts, there has been u rapidly expanding 
literature directed towards characterizing the u.v. 
spectral irradiance at the ground (Bener, 1972; Green 
et ai, 1974; Shettle and Green, 1974; Halpern et al., 
1974 ; Cutchis, 1974 ; Dave and Halpern, 1976 ; see also 
additional references listed in these works). As a result 
of these and subsequent works the problem of estimat- 
ing the solar spectral irradiance at the ground under 
clear sky conditions is fairly well advanced. However, 
the problem of allowing for the influence of clouds 
upon the u.v. spectral irradiance has thus far received 
only limited attention (Halpern et al., 1974 and Nack 
and Green, 1974) and much work remains to be done. 

In this paper we will examine the relative effect of 
cloud layers on received u.v. and net solar radiation 
through use of radiative transfer calculations applied 
in the case of uniform, homogeneous cloud layers. The 
result of variation of cloud thickness, cloud height, 
droplet size distribution, surface albedo and the at- 
mospheric ozone thickness upon u.v. and total solar 
radiation will be examined. The results presented will 
provide insight as to the relative influence of cloud 
cover on the u.v. and the net solar radiation. Our hope 
is that this insight will enable us to utilize the much 


* Supported in part by NASA Contract NAS9-15114. 


larger body of information available on net solar 
radiation reaching the ground in the presence of clouds 
to predict the u.v. radiation reaching the ground for 
which only limited data is available. 

2. COMPUTATIONAL PROCEDURE 

Received radiation at the bottom of the atmosphere 
is both direct and diffuse. Whereas, the calculation of 
the direct component is a simple procedure, calcu- 
lation of diffuse atmospheric radiation requires in- 
volved radiative transfer theo.y. A variety of radiative 
transfer methods have been applied for calculation of 
diffuse radiation in clear and cloudy atmospheres 
(Hunt, 1971). Both approximate calculation pro- 
cedures and more exact methods involving much 
computational effort are possible. The discrete or- 
dinate method for radiative transfer is an approximate 
solution, introduced originally by Chandrasekhar 
(1950), which gives good computational accuracy of 
fluxes with relatively rr, inor computational labor. The 
discrete ordinate radiative transfer procedure has been 
applied in recent years to cloudy and hazy atmos- 
pheres (Yamamoto et ai, 1971; Nack and Green, 
1974 ; Liou, 1976), and an adaptation of the discrete 
ordinate calculation described by Shettle and Green 
(1974) will be used in this study. 

The solution is applied to a vertically inhomo- 
geneous model atmosphere by dividing the atmos- 
phere into horizontally homogeneous layers and 
applying appropriate boundary conditions between 
the solution for various layers. For each layer, the 
transfer properties are determined by the layer optical 
thickness, the single scattering albedo and the scatter- 
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ing phase function, all of which may be calculated 
through knowledge of the molecular and aerosol 
scattering and absorption within the layer. 

Application of the discrete ordinate method as used 
in this study for calculation of u.v. spectral irradiance 
has been previously described by Shettle and Green 
(1974). Spectral irradiance values for a 2nm band 
interval are determined at specific wavelengths with 
absorption coefficients for ozone appropriate to this 
band interval being used (Green, 1966). The influence 
of the atmosphere on the u.v. spectral irradiance may 
be defined in terms of a spectral flux transmission 

where Fi„(A|,/<o) >s the downward spectral irradiance 
received on a horizontal surface at the bottom of the 
atmosphere, f(ki) the extraterrestrial solar spectral 
irradiance and Hq the solar zenith angle. 

Determination of the total received solar energy 
requires integration over ail solar wavelengths reach- 
ing the ground. By using broad spectral regions for the 
wavelength integration, computational effort can be 
decreased. Band weighted absorption parameters ap- 
propriate to the given spectral intervals must then be 
used. Fine wavelength integration intervals would be 
expected to produce more accurate results. However, 
for application of a more approximate radiative 
transfer solution such as the discrete ordinate method, 
the computation is intrinsically imperfect and the 
improvement of results by use of finer wavelength 
integration is limited. For this study the solar spectrum 
was divided so that the 0.7, 0.82, 0.94, 1.1, 1.38, 1.87 and 
2.7 nm water vapor absorption bands were included as 
separate spectral regions. The region from 0.34 fim to 
0.7 //m was further divided into four spectral regions of 
equal incident solar energy. For each spectral region, 
appropriate coefficients for molecular and aerosol 
scattering and gaseous absorption were determined. 
The coefficients for molecular and aerosol scattering in 
each spectral region were chosen to give values of 
scattering thickness equivalent to weighting by the 
spectral intensity of incident solar radiation. Absorp- 
tion coefficients for ozone, water vapor, CO 2 and 
other gases in the visible region were computed from 
the data given by Selby and McClatchey (1975). 

The flux transmission for the integrated solar energy 
may be given as 


F/s(/io) — 


S,-F,i(AA,) 

A'o'S 


( 2 ) 


where F,i(AA,) is the received downward energy in 
each spectral band for which a radiative transfer 
calculation is performed and S is the solar constant. In 
addition to the ten visible spectral intervals given 
above, the irradiance in the sub-340 nm region was 
determined by integrating over the spectral irradiance 
that was separately calculated for six u.v. wavelengths. 

Since the integrated solar irradiance is being de- 



Fig. 1. Comparison of the received diffuse solar flux calcu- 
lated by the wide wavelength interval, discrete ordinate 
method used in this study with the more exact calculations by 
Dave and Braslau (1975). Results for their Cl and DI-ST 
models are shown. Model Cl is a cloud free atmosphere with 
nonabsorbing particulates, and for Model DI-ST there is a 
stratus cloud layer and absorbing particulates. 


termined both with use of an approximate radiative 
transfer technique and with broad spectral integration, 
a comparison between the diffuse irradiance as calcu- 
lated by the method of this study and the more exact 
solution of Dave and Braslau (1975) was undertaken. 
The results are given in Fig. 1. For the calculations 
shown, the model atmosphere used for the discrete 
ordinate solution was set up to be the same as in Dave 
and Braslau’s (1975) Cl and Dl-ST models. The 
discrete ordinate calculation was applied with an eight 
layer atmospheric model. In the case of the cloud free 
Cl model, the values for diffuse transmission of solar 
energy through the atmosphere agreed for all zenith 
angles. For the DI-ST model with a 3.57 optical 
thickness (0.55 fj) stratus cloud, the discrete ordinate 
calculation produced values of received diffuse solar 
energy that were a few per cent low for small zenith 
angles and a few per cent large for the extreme zenith 
angles. 


3. MODEL 

For all calculations in this study, aerosol scattering 
parameters were determined using the “Haze L” size 
distribution by Deirmendjian (1969). Unless otherwise 
noted, an aerosol refraction index of 1. 50-0.01 f and a 
total particulate optical depth ofO.125 at 590 nm were 
used. The particulate height distribution was that 
given by Shettle and Green (1974). 

The water vapor and mixed gas profiles for the 
model atmosphere is the representation given by 
McClatchey et al. (1970) for midlatitude summer 
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conditions. The ozone vertical profile given by Shettle 
and Green (1974) with a total ozor . thickness of 
0.318 cm is used unless otherwise noted. 

Cloud layers will be defined by the total optical 
thickness at 590 nm and by the height for the top and 
bottom of the cloud layer. The phase function for cloud 
scattering was computed using the Cl cloud size 
distribution ofDeirmendjian with a refraction index of 
1.34 +0.0i. Legendre polynomial coefficients for use in 
the discrete ordinate transfer solution were calculated 
directly from the cloud scattering phase function. 
Unless otherwise defined, a cloud with a base at 3 km 
and top at 4 km will be used for model calculations. 

4. RESULTS 

The relative influence of cloud cover on received u.v. 
and total solar energy will be presented in terms of the 
ratio of u.v. and total solar energy flux transmission 
through the atmosphere. 

= < 3 ) 

Values of Ri were calculated as a function of cloud 
thickness, cloud height and surface albedo. Results for 
five u.v. wavelengths are to be presented 2<J0, 300, 
310, 320 and 330 nm. 

The behavior of the u.v. to total energy transmission 
ratio as a function of solar zenith angle in the case of 
zero surface albedo is shown in Fig. 2. The decrease of 
u.v. spectral irradiance as the zenith angle increases is 
primarily due to ozone absorption. The actual values 
of the u.v. total energy transmission ratio will be 
strongly dependent on the atmospheric ozone thick- 
ness, and the values given in Fig. 2 apply only for a 
thickness of 0.318 cm of ozone. 

As seen in Fig. 2, the presence of a cloud layer of 
increasing thickness acts to increase the ratio of the u.v. 
to solar energy flux transmission through the atmos- 
phere. One possible factor is that both the cloud 
optical thickness decreases and the droplet scattering 
phase function becomes more sharply forward-peaked 
for shorter wavelengths. Values of cloud optical thick- 
ness, Tf, and the asymmetry factor of the cloud 
scattering phase function, g, are given in Table 1 as a 
function of wavelength The transmission through an 
isolated cloud layer is thus greater for shorter wave- 
lengths. However, it was found that changes in the 
cloud droplet distribution did not appreciably effect 
the ratio of u.v. to solar energy flux transmission 
through a cloudy atmosphere. Even when the values 
given in Fig. 2 are recomputed using a constant cloud 



l ig. 2 The ratio of total transmission, diffuse plus direct, for 
u.v, radiation and total solar energy as a function of solar 
zenith angle. Values for three different cloud thicknesses, at 
each of live wavel.;ngths, 290, 300, 310, 320 and 330 nm are 
shown. The scale o.i the left applies for all wavelengths except 
290 nm for which the scale on the right is used. In this case a 
surface albedo of 0.0 was assumed. The values shown are for 
an ozone thickness of 0.318 cm. 


optical thickness and scattering phase function with 
wavelength, the results are not significantly changed. 

The behavior of the u.v. to solar energy flux 
transmission ratio with increasing cloud thickness is 
due primarily to the interaction of the cloud layer 
scattering with that of the surrounding atmosphere. 
The scattering thickness of the cloud-free atmosphere 
is much greater at u.v. wavelengths than for visible 
wavelengths. In the case of zero surface albedo as in 
Fig. 2, there will be significant reflectivity by molecular 
scattering beneat h the cloud layer for the shortwave 
radiation but very little reflectivity under the cloud 
when considering the integrated solar radiation. A 
doubling interaction between the cloud layer re- 
flectivity and scattering by the underlying atmosphere 
results in a smaller decrease of the surface u.v. ir- 
radiance with increasing cloud thickness. The effect is 
moderated by the u.v. radiation incident on the cloud 
top being more diffuse than the integrated solar 
radiation, but not sufficiently moderated to give a 
relatively greater decrease of the surface u.v. irradiance 
with increasing cloud thickness. 

It can also be seen from Fig. 2 that the behavior of 


Table 1. Cloud optical thickness t, and scattering phase function asymmetry factor 9 as a function of wavelength 
2 290 310 330 405 490 590 690 790 920 1110 1400 1900 3050 


9 


4.56 4.69 4.78 4.91 4.96 5.00 5.03 5.07 5 11 5.16 5.25 5.38 5.79 

0,860 0.860 0.859 0.858 0.856 0.853 0.850 0.846 0.841 0.834 0.823 0.802 0.768 
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Realistic cloud cover, however, is not uniform and 
homogeneous. For broken cloud cover, the ratio of the 
instantaneous received u.v. and net solar radiation will 
change rapidly. Since the u.v. radiation incident on a 
cloud layer is much more diffuse than the visible 
radiation, a partial coverage cloud which blocks the 
direct solar radiation will decrease the total received 
visible flux to a greater extent than the total received 
U.V. flux. Over time periods sufficiently long to average 
such fluctuations, however, it might be expected that 
the relative decrease of the received u.v. and net solar 
radiation would be similar to that predicted for 
uniform cloud cover of an appropriate thickness. 
Either much more elaborate radiative transfer calcu- 
lations or a body of experimental data would be 
required to determine if this is the case. 
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ABSTRACT 

The purpose of this study is to quantify for the 
needs of photobiology the influence of clouds upon the 
ultraviolet spectral irrsdiance reaching the ground. 
Towards this end we firs, develop analytic formulas 
which approximately chr^iacterize the influence of clouds 
upon total solar radiation. These may be used in con- 
junction with a solar pyranometer to assign an effective 
visual optical depth for the cloud cover. We also de- 
velop a formula which characterizes the influence of 
the optical depth of clouds upon the UV spectral irra- 
diance in the 280-340 nm region. Thus using total 
solar energy observations to assign cloud optical prop- 
erties we can calculate the UV spectral irradiance at 
the ground in the presence of these clouds. As inci- 
dental by-products of this effort, we have found con- 
vient formulas for the direct and diffuse components 
of total solar energy. 


1. INTRODUCTION 

There is now great concern that anthropogenic modif icat .ons of the ozone 
layer might alter the ultraviolet (UV) radiation reaching the ground. This 
concern has led to a rapid expansion of the literature on the biological im- 
pacts of UV radiation and on the quantitative aspects of the UV spectral irra- 
diance reaching the ground [Bener (1972), Green, Sawada and Shettle (1974) 

(GSS), Shettle and Green (1974), Halpern et a l., (1974), Dave and Halpern 
(1976)1. Unfortunately clouds are a large source of uncertainty in such wor)c. 
Naclt and Green (1974), Halpern, Dave and Braslau (1974) , and Spinhirne and 

Green (2978) have carried out radiative transfer calculations of the influence 
of cloud cover on the received UV irradiance for wavelengths in the 280-340 nm 
range, at various solar zenith angles, and various ground reflectivities. To 
facilitate the practical use of such calculations in photobiological applica- 
tions, we attempt, in the present work, to express the results of these earlier 
works and our supplementary calculations in terms of convenient analytic forms. 

The strategy of the present effort is to use observations of total solar 
energy to infer the optical properties of clouds. Thus the reduction of total 
solar irradiance from the expected cloud-free levels can be used to estimate 
the effective cloud optical depth, the critical parameter needed to improve UV 
calculations. Characterizations of total solar energy will be discussed in 
Section 2. 


2. CLEAR SKY TOTAL IRRADIANCES 

For reference in our cloudy sky work, it is valuable to have approximate 
representations of the clear sky irradiances for various sun angles and ground 
reflectivities. Using Beer's law the monochromatic downward component of the 
transmitted (t) solar irradiance may be calculated using 

= pH(A)exp - w^^ (y) k^ ( X ) (1) 
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where HU) is the extraterrestrial solar spectral irradianoe, 9 is the soler 
zenith angle, i « 1, 2, 3 denote air, particulates and other absorbing compon- 
ents, Wi(y) denote the altitude dependent vertical thickness and k. (X) the wave- 
length dependent attenuating coetCicient. To allow approximately ror the fact 
that the earth is round, we will use a genera titration of m » cos9 given by 

Mi « i'(M^ + K^)/(l + (2) 

where Kj are species dependent dimensionless parameters. The equivalent gen- 
eralization of sec9 » p”i is o =* seq9 - (GSS, 1974) . 

To calculate t)'e "transmission" for the downward component of the total 
transmitted (t) solar imdiance, we use 

-E.w.k./M. 
li S AA^H(Xj)e ^ i ^ ^ 

'^Tt “ ^ (Ts 

where S = E AX.H(A.). The wavelength intervals used in the present calcula- 
tions are given in^Spinhirne and Green (1978). 

We may represent the results of our total solar transmission calculations 
to a good approximation by the "random" model which Mayer (1947) and Goody 
(1952) developed for infrared band transmission. Thus we may use 


T 


Tt 


exp 


gg 

(l+So)’’ 


(4) 


where the k parameter used in a has been determined by the rule 

^ = ‘Vr + ^ V 

Where k- * 1,8x10 , k * 0,3x10 ana and t are the Rayleigh and particu- 

late optical depths at"550 nm. Figure l^shows "the calculated values (points) 
using Eq, (3) and the analytic representations (solid lines) using Eq. (4) for 
the various aerosol optical depths with parameters given by 


a 


“1 + 


“2"p 


( 6 ) 


with a, = 0.2138 a, = 0.8848 and 3 = 0.1717. It is seen that the random 
model transmission formula with this assignment of parameters works quite well. 
The result is not entirely unexpected if one considers the physical analogies 
involved in the calculation of infrared band model transmission and the calcula- 
tion of total solar energy transmission (see Green and Wyatt, 1965) . 


To calculate the downward diffuse (d) irradianoe (skylightl) we adapted 
the discrete ordinate radiative transfer code of Shettle and Green (1974) to 
obtain the "diffuse transmission" for total solar radiation given by 


T 


Td 


A AA.F-KA.) 
__ i- 


(7) 


We have represented the calculated results analytically using 




K exp - 


yq 

(l+6a)*^ 


( 8 ) 


with the parameters given by 


and 


K = + 


'^2'^p 


Y = <5 = Yi + Yz^p 


(9) 
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(10) 


with K, » 0.0716 K- » 0.8642 y, » 0.283 y~ » 1.302 . The lower c <rv»8 in 
Fig. l'‘’«how the fits to the numerical calculations (points) . The fits arc sur- 
prisingly good considering the fact that the physical processes are so dxrfer- 
ent from those implied in the random model. The fits can be adjusted to well 
within 1% if we relaxed the specification that y “ 6. 


To allow for the influence of ground albedo (liambertian reflectivity) we 
may use a device of Shettle and Green (Eq.43). T)ii8 is based upon the assign- 
ment of a reflectivity parameter r which is the fraction of the total radiation 
reflected by the surface which is scattered back to the surface by the atmo- 
sphi ;e. Thus for any albedo the global transmission for solar radiation may be 
represented approximately by 


'Tg 


^Tt ^Td 

aT" 


( 11 ) 


where r * 0.124 and and given by Eqs. (6) and (8) respectively. 

3. CLOUD EFFECTS UPON TOTAL SOLAR RADIATION 


Dave and Braslau (1975), Liou (1976), and Spinhirne and Green (1978) have 
calculated the total solar radiation penetrating various cloud thicknesses. 

Here we consider the ratio of the downward component of solar radiation in the 
presence of a homogeneous cloud layer of optical thickness t (at 550 nm) to the 
total downward solar radiation for a cloudless sky. Thus we characterize the 
transmissivity of the cloud layer by the ratio 


R ( I , n ) 


r .. f^i(aa^,h,t) 

T/. F^(AA^,U,0) 


( 12 ) 


where F.4 -(AX.,m,t) is the total received downward energy in each spectral band 
used in^the Jadiative transfer calculation. 


We have obtained good fits to the results of our radiative transfer calcu- 
lations with the surprisingly simple formula 


R,j,(ifM.A) 


T (1-rj^A) 

exp 2 2 >7 

[a‘^(Vi-lio) + ct(l-r2A)l^ 


(13) 


where a = 16.0, y = 0.240, c = 9.13, r, = 0.922 and r, = 0.952. Figure 2 
illustrates the fits of vs. y for t = 2, 10, and 30 when A = 0.1 and A = 0.5. 

We see that Eq. (13) provides an analytic representation of the three di- 
mensional hypersurface R^(i,y,A) which expresses the variation of the total 
solar radiation ' h sun angle (through y) ground albedo (A) and cloud thick- 
ness (t) in ter: *jt only five adjusted parameters. Let us next consider the 

corresponding representations for spectral irradiance in the 280 to 340 nm 
region. 


4. ULTRAVIOLET SPECTRAL IRRADIANCES FOR CLEAR SKY 


The ultraviolet spectral irradiance problem has the wavelength (A) and the 
ozone optical depth (13) as important additional degrees of freedom. The down- 
ward transmitted spectral irradiance at the ground may be computed simply using 


B^.(A,0) = uH(A) exp 



(14) 


where H(A) is the extraterrestrial solar irradiance, t,, t_, and t, denote the 
wavelength dependent Rayleigh, aerosol and ozone optical depths, respectively, 
and y, yj^, y^ and y3 denote the cosf) and generalized cosine functions for these 
species. If we now specify the optical depths of the three components at a 
standard wavelength (Aq = 300 nm) and their dependence on wavelength, we can 
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■imply cAlculate the downward transmitted component of the sun's UV irradiance. 
The apectral dependence of the optical depths for ozone, air and a representa- 
tive aerosol model in the wavelength region of interest may be represented 
analytically by 






and 


r r . g ^ ■ 1 

30 + exp{A-X^)/d _ 


where K * 300 nm, “ 1.221, v = 4.27, t 
d “ 8 nm, and 730 >■ “3.2, Of co 8 rse, 
conditions. 


Tp an 


(15) 


p - 0.365, Vp » 0.58, 3 = 0.0439, 
d Tqz are dependent upon atmospheric 


Skylight or the diffuse irradiance is particularly important in the ultra- 
violet. Dave and Furukawa (1966) have calculated UV skylight for the case of 
a purely molecular atmosphere. To deal with the turbid atmosphere. Green, 
Sawada and Shettle (GSS, 1974) developed a semi-empirical analytic representa- 
tion to fit the synthesis of 10 years of UV spectral irradiance measurements 
by Bener (1972). In addition, Shettle and Green (1974) have used a fast and 
inexpensive radiative transfer calculational technique to obtain the UV dif- 
fuse component in a turbid atmosphere. Dave and Braslau (1975) have also com- 
ple»-«(’ UV radiative transfer calculations using a much more refined computation 
whic I ,an serve as a standard of accuracy in atmospheric radiative transfer 
calculations. Ihese sets of calculations can also be approximately represented 
by the GSS empirical formulas with minor parameter adjustments,, 


The GSS semi-empirical formula for the diffuse spectral irradiance has 
appeared in the literature with several minor variations and with somewhat dif- 
fering parameter sets. It is believed that their analytic representations of 
radiative transfer calculations or experimental data can be improved from the 
25% level of accuracy to the 5% level by using molecular band models (Green 
and Wyatt, 1965) rather than Beer's law as a point of departure. However, 
pending the completion of this work, we recommend using the version and para- 
meters of the GSS formula in Green, Mo and Miller (1974) . To correct for 
ground reflectivity upon clear sky global spectral irradiance, we use the 
ground-air ground-reflectivity function of Shettle and Green (Ea, 46) . 

5. CLOUD EFFECTS UPON UV SPECTRAL IRRADIANCES 


To represent cloud effects upon UV spectral irradiances we use the same 
basic methodology as that useu to represent cloud effects on total solar radi- 
ation. Thus we assume that the ratio of the global cloudy sky spectral irra- 
diance to the global clear sky spectral irradiance (for A = 0) may be repre- 
sented by a slight generalization of the form 


^uvg 


exp 




{a^[(M-MQ)^ + + cr}^ 


(16) 


where a = 3 q[Tj^/(Tj^ + and a^ = 25.14, = 0.4181, k = 0.732, c = 16.26. 

Figure 3 illustrates the fits of Ruvg t = 2, 10 and 30 and for 

several wavelengths. The magnitudes of the results are quite good although 
the small variations with p are not precisely described. However, the simpli- 
fied radiative transfer calculations themselves m, ’it not be accurate on such 
fine points. 


We have explored the dependence of the R,ivg upon j round reflectivity. We 
have had reasonable success with the simple device of replacing x in Eq. (16) 
by xd-rA) with r = 0.69 (see Eq. (13)1. Undoubtedly we can do better by let- 
ting r be wavelength dependent but since UV ground reflectivities are low (Fu- 
rukawa and Heath, 1973) this simple correction should handle many situations. 
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APPLICATIONS 


The basic strategy of this applied study is to take advantage of the large 
body of solar radiation measureitents to assess the optical property of clouds. 
Thus the local user would use observations of ground level total solar irradi- 
ances together with Tto* given by Eg, (11), and R^, given by Eg, (13), to estab- 
lish the c values of clouds in the visible region. Here we can take advantage 
of the fact that Eg, (13) is algebraically invertible. First we cast Eg. (13) 
to the random model form 

R exp - (i/[j^)/ll + ( 1 / 12 )]*^ • (17) 

Then by algebra we find the visible optical depth, i, using 



Next we translate visible ground albedo into UV ground albedo using the data of 
Heath and Furukawa U973). Then using the CSS formula for the clear sky, to- 
gether with Eg. (16), we calculate the UV spectral irradiance at the ground 
from a cloudy sky. Wo hope in time to relate our results to spectral irradi- 
ance observations from satellites (Heath, Matecr and Krueger, 1973). 

It should be noted that at ■ 300 nm the cloud transmission ratio [Eg, (16)1 

is sensitive to the ozone concentration in the cloud (Spinhirne and Green, 1978). 
The characterization of this sensitivity remains to be established. 

When the cloud thicknesses are not uniform we must first establish a dis- 
tribution function for the optical depth. This can be done by using a time 
trace of the total solar radiation to infer [via Eg. (13)] the time variation of 
! over the observing site. Then using those 1 in conjunction with Eg. (16) 
and GSS formula for the clear sky UV spectral irradiance we may determine a 
distribution function for the UV spectral irradiance for a cloudy sky. 

We are currently involved in practical tests of the foregoing procedures 
by independent measurements of total solar irradiances and UV spectral irradi- 
ances when clouds are present. From our preliminary investigations wo are con- 
fident that the procedure should work reasonably well whoa the horizontal scale 
of inhomogeneities is small compared to the vortical scale. It is not yet clear 
how the method would work for the case of widely scattered clouds but such 
studies arc also underway. 


7. SUMMARY AND CONCLUSIONS 

We have defined a procedure for translating observations of total solar 
irradiances when clouds are present into estimates of UV spectral irradiances. 
The procedure can bo implemented with the aid of a modern desk top computer. 

The essential aspect of our procedure is the use of semi-empirical representa- 
tions of diffuse radiative transfer obtained by outputs of complex radiative 
transfer calculations made with largo computers. Our work indicates that these 
outputs lend themselves to analytic representations which arc far more conveni- 
ent in application than the tabular outputs themselves. For example, once the 
analytic representations of the UV spectral irradiance vs. a, etc. are 
available, one can by the methods described by GMM calculate daily UV dose, 
monthly dose, annual dose, etc. One can also integrate over various measured 
or assumed spectral I'esponse function (action spectra) . The present work is a 
feasibility study which indicates that it is possible to implement such metho- 
dology. 
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ANALYTIC MODEL APPROACH TO THE INVERSION 
OP SCATTERING DATA 


Alex B. S. Green and Kennath F. Klenk 
University of F lor ids 


Me apply an analytic model approach which has _ been 
developed in nuclear studies to several simple atmospheric 
inversion problems. We illustrate by past work on the 
solar aureole that this method gives a sharp determination 
of aex:>sol size distribution parameters. We show that this 
analytic approach, together with ground level point sampling 
data measurements, may be used to infer information on the 
tropospheric ozone protiie. 


I. ATMOSPHERIC AND NUCLEAR OPTICS 

Many of us who are now involved in atmospheric inversion 
problems were previously involved in analogous problems in other 
disciplines. As is natural, we try to bring to bear the experience, 
sense of aesthetics, or prejudices, if you will, which we have 
acquired in these other fields. The beauty of this conference as 
it is developing following some of the earlier papers is the sense 
of open-mindedness which has emerged. It is as if this conference 
has said, "Let a thousand flowers blossom." 

In my own (Green) case, my main prior involvement with inver- 
sion problems has been in connection with two nuclear physics endeav- 
ors based largely upon scattering data — {%) inferring the nature of 
the fundamental interaction between neutrons and protons, and 
(2) inferring the detailed nature of the nuclear potential mani- 
fest in the shell and optical models of the nucleus. Let me use 
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ttM first. problm as an illustration of hov undsrstandin? is 
advanced in scattering inversion problems. 

When the neutron was discovered in 1932 > the fundamental prob- 
lesi in nuclear physics became that of inferring the basic force 
between neutrons and protons. The approach followed has bean to 
perform scattering experiments, l.e., fire neutrons or protons on 
hydrogen targets, and examine the emerging angular distributions 
(phase functions) and polarizations at various energies (wave- 
lengths) of the outgoing particles. The hope was to be able to 
test various proposed two-lsody potentials which when inserted into 
the Schrddinger equation or Dirac equation might account for these 
data within statistical error. This was the main line of approach 
in nuclear physics until the early 1960s when the only phenomeno- 
logical models which could fitf the 0 to 400 MeV array of scattering 
data and auxiliary data such as the properties of the bound two- 
body system (the deuteron) were exceedingly complex, requiring 
as many as 40 adjustable parameters in their description. 

h brea)ct)irough came in the mid*1960s when the discovery of the 
u, p and n mesons by particle physicists led to the revival of 
meson theory of nuclear forces initiated by Yukawa in 1935. With 
the additional physical constraints of meson theory, it suddenly 
became possible to fit the scattering data with one boson exchange 
model requiring only five to ten adjustable parameters, rather 
than the 40 parameters of purely phenomenological models. Although 
the final story is not yet told, the nuclear physics community, 
since 1967 (Refs. 1 and 2), has felt a great aesthetic sense of 
relief that the fundamental law of nuclear physics is not as mon- 
strous as it had appeared to be in the (sarly sixties. 

Thus, as some of the earlier speakers have already suggested, 
it is the additional physics, physical judgment, and physical infor- 
mation which one brings to bear with the scattering data which will 
often determine the success and utility of an Inversion scheme. 
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Studies leading to the nuclear shell and optical Models 
(Ref t 3) are even More analogous to the atMOspheric inversion prob- 
ISM. For over 40 years, experlMents have been perfomed in irttich 
various nuclear particles accelerated froM 1 HeV to 100 GeV energy 
range are scattered from various nuclear targets. Many lead to 
optical -type angular distributions (phase functions ) , polarisations, 
scattering, and absorption cross sections. Many of these data 
patterns can be accounted for by assuming a complex energy depen- 
dent nuclear potential (complex wavelength dependent index of 
refraction) . Even the terminology of this subject, such as "the 
cloudy crystal ball model," reflects the light scattering analogy. 
Mow nuclear opticians, like atmospheric opticians, divide up into 
a sc)tool concerned with average gross properties and a school con- 
cerned with statistical fluctuations*. Both groups have greatly 
enriclied the subject, although, as in the light scattering case, 
the communications between the schools has not always been the 
best. 


My own specialized pursuits of atmospheric optics (apart 
frem a stint in World War II) began in 1959 just after an intensive 
involvement with nuclear optics (Ref. 4). In these pursuits, I 
have mostly used the gross structure-nuclear optical modelers 
approach. The style here has been to use analytic models whose 
parameters are determined by nonlinear least square adjustment to 
experimental data. Then we look at the systematics of the param- 
eters with the ultimate objective of relating them to more funda- 
mental physical parameters, e.g., those in the basic nuclear force. 
X would like now to illustrate this nuclear optical approach with 
a few simple-minded attacks on some atmospheric optics problems. 

XI. AUREOLE STUDIES 

Deirmendjian and Sekera (Refs. 5 Md 6) very early recognized 
the Isportance of Mie pa.rticles in the theory of the solar aureole. 
Their work was motivated by an attempt to account for some reported 
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•noMlously high transmission of the ultraviolet part of direct 
sunlight. It is interesting to note that D. S. Saxon* who col- 
laborated with Deitmendjian and Sahara on scattering from dielectric 
spheres (Ref. 7) during the same time frame* played an Isqportant 
role in the development of the nuclear optical model. 

Our %«ork* on determining aerosol sixe distributions from solar 
aureole Intensities* has used a type <~>f atmospheric modeling and 
single scattering theory which we first used in a satellite ozone 
sounding analysis (Ref. 8). That the skylight in the neighborhood 
of the sun can be used to good advantage may be surmised from the 
optical theorem of nuclear and atomic scattering theory and 
bistatic radar analysis (Ref. 9) . These works indicate jthat 
whereas backscatter cross sections (180° scatter) vary in a 
complex manner with particle shape 'and index of refraction* for- 
tiard scatter cross sections are primarily determined by the volume 
of the particle. This property was first utilized in the doctoral 
theses of Adarsh Deepak (Ref. 10) and Barton J. Lipofsky (Ref. 11) 
which* among other things* involved photographic and photoelectric 
studies of the solar aureole in the visible region (Ref. 12). More 
recently* these studies have been extended into the ultraviolet 
(Refs. 13 and 14) . 

Or. Deepak has already described some of the features of this 
work in his talk on the Aureole Isophote Method and in connection 
with his Stratospheric Aerosol Photographic Experiment (SAFE) 
proposal. Let me add some words here on the advantages of this 
type of experiment. 

In Fig. 1* we show the densltometric traces of a measured 
photographic aureole made with the obscuring disc technique. The 
traces labeled r emd 1 are through the solar almucantar whereas 
the traces labeled b' and t are below and above the sun in the sole.r 
meridian. The beauty of this experiment is the reference intenslcy 
provided by the solar disc, so in essence we have on our measuring 
medium a comparison between diffuse sky intensity and the direct 
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Fig. 2*. Bxamp2es of oversize distribution. 

Fig. 2b. Corresponding exampies of distribution. 

solar intensity, attenuated by lO** by a Neutral Density 4 (ND4) 
filter. Thus far, in our work at the university of Florida, we 
have only exploited the information content in the almucantar 
trace of the solar aureole. 

For' the convenience of analysis, we use an analytic size 
distribution character! sed by the cumulative distribution function 
(see Fig. 2) . 
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N* hav* calculated a library of normalizad phase function for 
this two-parameter size distribution function. Figure 3 illus- 
trates several examples. Figure 4 illustrates some recent work 
at 640 nm showing the sharpness in the determination of V obtained 
by this type of analysis (Ref. 14). 

fipectral turbidity measurements, such as determined by a 
multi-channel Sun photometer, also give information about the 
particle size distribution (Ref. 15). Green and Sawada (Ref. 16) 
have determined the relative spectral turbidities associated with 



(b) V » 5. 

Fig. 3. Concluded. 
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Fig. 4. Measured aureole intensicy* along the almucantar 
compared with calculations for three v variations in the aerosol 
sise distribution. (From Ref. 14.) 

aerosol eize distributions characterized by Eqs. (1) and (2) . 

Thus* aureole data augmented by spectral turbidity data can lead 
to a very sharp determination of v, as well as N^. 

We have considered the problem of multiple scattering in the 
solar aureole (Ref. 17 ) , particularly in the ultraviolet where 
Rayleigh scattering becomes so important. Figure 5 shows an inter- 
comparison of three calculations (Ref. 18) : (1) a single scat- 

tering treatmentf (2) a multiple scattering calculation using Monte 
Carlo techniques, and (3) a multi-channel calculation. The cal- 
culation indicates that multiple scattering can be quite signifi- 
cant, particularly in the ultraviolet. 

Zn subsequent woric, McPeters emd Green (Ref.. 14) have found 
that Rayleigh scattering is the source of most of this multiple 
scattering in the solar aureole except at high particulate optical 
depths. Accordingly, they have proposed an analysis technique 
which uses single aerosol and Rayleigh scattering augmented by 
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Fig. 5. Absolute intensity as a function of detector zenitli 
angle for two different Hie scattering functions as calculated by 
three separate methods. The solid lines represent the results of 
th« multi-channel calculations, the dots are the Monte Carlo 
results, and the dashed curves are the single scattering calculations 
Mote the break in the scale to avoid superposition of the two sets 
of results. (From Ref . l£.) 
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iNiltipl* Raylalgh scattering as determined by the use of the tabJ.es 
of Coulsoni Dave«and Sekera (CDS) (Ref. 19). Figure 6 is an illus- 
tration of a fit to data based upon such an analysis. 

XXX. DXFFUSE TO DXRECT RATIO MEASUREMENTS 

Xn connection with our ultraviolet studies in support of the 
Clisiatic Xsg>act Assessment Program (Ref. 20), we have attempted to 
realistically characterize the diffuse solar radiation or sky 
radiation in the ultraviolet. As it turns out, sky radiation in 
the ultraviolet is often of greater biological consequence than 
direct sunlight. Green, Sawada and Shettle (Ref. 21) have devel- 
oped an approximate analytic formula which describes the diffuse 
spectral irradiance in the ultraviolet region by adapting a single 
scattering analysis to the systematics of Bener's experiments 
(Rsf. 22) and to theoretical calculations of Shettle and Green 
(Ref. 23). 



Fig. 6. Normalized ultraviolet aureole intensities along the 
almucantar compared with modified single scattering calculation, 
with and without the CDS Rayleigh multiple scattering correction 
factor, and compared with a Monte Carlo full multiple scattering 
calculation. (From Ref. 14.) 
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Hor« r«c«ntly» Chal and Graen (Raf. 24) hava raeognifod tha 
■wrlta o£ maasuremanta of tha ratio of diffusa to diract apactral 
irradianca aa a simpla indicator of atmoapharic optical propartxo. 
Tha ratio method is analogous in soma raspacts to our aureola 
method except that there wa must ajttenuate tha diract solar inten* 
sity by four orders of magnitude to compare it to tha sky intensity. 
Figure 7 is an example of measuramants of the total sky irradianca « 
t)ia diract irradianca/ and the ratio. The diffusa and direct both 
vary vary markedly, and reflect the fluctuations in the extra- 
terrestrial solar spectral irradianca and in the ozone extinction 
coefficients as well as in the wavelength dependence of instrumental 
sensitivity. However, tha diffuse to direct ratio is only slowly 
varying but still sensitively depends upon such interesting charac- 
teristics as atmospheric particulat^ loading, ground albeda and 
sky cover. This ratio method avoids problems associated with the 
difficulty of absolute spectral irradiance measurements which, at 
this time, are limited to about 8% in the ultraviolet region. 

It should be remarked that Herman, et al. (Ref. 25) have 
theoretically examined such a ratio method in the visible region 
in an attempt to estimate the imaginary part of the index of 
refraction of atmospheric aerosols. 

He shall next consider another example of inverting of optical 
data with the aid of the analytic modeling approach and auxiliary 
Information, such as may be obtained with simple ground-based 
instruments. 

IV. THE INVERSION OF THE LOW ALTITUDE OZONE PROFILE 

Tropospheric ozc> :< is a constantly varying atmospheric com- 
ponent which changes with the season of the year, location, and 
time of day. Both photochemical and stratospheric transport pro- 
cesses are important sources of tropospheric ozone and are respon- 
sible for the strong diurnal, seasonal, and spatial variations. 
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fig. 7. Total, diffuse, direct, and the ratio of diffuse to 
direct components of solar irradiances as a function of wavelength 
as recorded on 31 July 1975 at around 13:50 EST. The magnitudes 
of the irradiances are in relative scale; they are proportional to 
the photomultiplier output in 100 pv. (From Ref, 24,) * 


A systematic program of balloon-borne ozonesonde observations 
has provided valuable data on the altitude structure of ozone which 
can be used as a guide in constructing a versatile analytic model 
(Ref. 26) . We consider in what follows the possibility of infer- 
ring the tropospheric ozone profile from diffuse to direct ratio 
measurements in the middle ultraviolet in conjunction with ozone 
point-sampling at the ground. We model the altitude profiles of 
the atmospheric components with analytic functions because such 
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a tachniqua slaipllflaB data invaraion and providaa a convanlant 
way of eonwunicating tha rasulting prof lias. 


Graan (Raf. 8) has nodalad tha stratospharle osona eolum 
danslty as a function of altltuda y with a distribution usad 
axtanslvaly In tha nuclaar studlas (tha so-eallad Wood-Saxon 
function) (Refs. 3 and 4) 


w(y) - 




'o/h) 


1 ♦ a‘y " 


(3) 


Bora w^ Is tha total ozone thlclcness and and h are parameters. 
Tha danslty profile p(y) Is given by 


P(y) 



(1 ♦ e 


^o/h 


La 


(i+*e<y- 


( 4 ) 


The parameter y^ Is the altitude at which the density function 

peaks and h scales the width of the distribution. Green (Ref. 2) 

shows how y^, h and w^ can be approximately Inferred from solar 

backscatter ultraviolet (UV) measurements.. In their recent 

analysis on ground level UV, Shettle and Green (Ref. 23) add an 

exponential term to this function to allow for the tropospheric 

ozone component. Here to characterize tropospheric ozone profiles 

which are concave, l.e., the density decreases with* altitude above 

the ground and then Increasing at the tropopause, we add a second 

term of the form of Eg. (3) or (4) with the parameters y ' and h' 

and where w_ Is now the sum of‘w + w '. 

T o o 

In Fig. 8, several profiles corresponding to various values 

of h' and p(0) are shown. Here y_* - 0 and w « 0.29 atm-cm> and 

o o 

y^ and h are set to 23 km and 4 km, respectively. Extreme con- 
cave profiles can be obtained as well as curves of almost constant 
density. Furthermore, convex profiles can be generated with two 
distributions by setting y^' to be a positive number. Convex pro- 
files are observed with greatest frequency In the summer months at 
latitudes above 4C°N. 
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Fig. 8. Ozone cionsity profiles for a number of different 
values of the parameter h^’ and ground density, with y^ « 23 km, 
h • 4 km and 1/^*0. 

‘Asauffling we know the ground level ozone density and also tltat 

the ozone profile can be represented by a two component model vflth 

y* • 0> we can test the sensitivity' of the diffuse to direct ratio 

to the parameter h'. In Fig. 9, we plot the ratio versus h* for 

p(0) “ 40f 70> 100 (yg/m^) for a wavelength of 300 nm. We take 

y. 23 km and h_ « 4 km. Furthermore, we have assumed the air 
o o 

and aerosol profiles of Shettle and Green (Ref. 23) which are based 
on Elterman's 1964 data (Ref. 27). The total aerosol optical 
depth is 0.411; the aerosol is characterized by the cumulative 
size distribution given by Eq. (1) with v •• 3 and a *0.03 uni 


the 
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Figr. 9. Diffuse to direct ratio dependence on the scale 
parameter h' for ground ozone densities of 40, 70, and 100 


aerosol index of refraction is l.S .Oli; a ground albedo of 0.05 
is assumed; and the Sun is directly overhead. The ratios are cal- 
culated by using the multiple scattering technique of Shuttle and 
Green <Ref. 23). 

The diffuse to direct ratio then can be used to infer a value 
of h* fr«n a set of curves as in Fig. 9 when used with point 
sampling measurements which give the ozone concentration at the 
ground. 

Other simple and inexpensive measurements can be made simul- 
taneously with the ratio and point-sampling measurements in order 
to further delimit the inversion and reduce the uncertainties. For 
example « aureole photography can be employed to infer the aerosol 
size distribution parameters (Refs. 10> 12, and 14) . Also, 
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■ulti-wav«l«ngth photometry provides valuable infomatlon on the 
aerosol optical thickness and its wavelength dependence (Refs. IS 
and 16) . 

The two most nebulous parameters are the aerosol single 
scatter albedo and the ground albedo. It is important to know how 
uncertainties in these parameters propagate through the inversion. 
The aerosol single scatter albedo will depend on the aerosol index 
of refraction and sixe distribution. By using aureole photography 
to pin down the size distribution, one can draw on experience from 
previous Investigations to fix the aerosol index of refraction 
within a range of confidence. For example, blstatic laser and 
aureole photography methods (Ref. 28) indicate that the index of 
refraction of a typical Gainesville aerosol is O.OOS ± 0.005 for 
the imaginary part and 1.50 i 0.05 for the real part. The diffuse 
to direct ratio is found to be relatively insensitive to this range 
of possible error. Similar considerations apply to the ground 
altMdo. In the theoretical model a ground albedo of 0.05 was 
assumed which is compatible with measurements by Furukawa and Heath 
(unpublished reports. 1973) of various natural surfaces for the wave- 
length region 310 to 380 nm. For example, they found that for scrub 
desert the ground albedo was 0.04 over the 310 to 380 nm region. 

For farmland. 70% tilled and 30% covered with vegetation., the ground 
allaedo was found to be 0.07 to 0.08 for the 310 to 340 nm region. 

Small errors in the ground albedo («0.02) do not significantly 
affect the calculations. Furthermore, once the ground albedo is 
known for a given location, the dally variations of h' can be 
determined, unless, of course, the surface c)ianges because of snow 
cover, cultivation, or the like. 

Zf the true optical depth is used in the inversion, then 
underestimating the ground albedo or the aerosol single-scatter 
albedo will lead to calculated ratios which are too small. 
Overestimation will lead to calculated ratios which are too large. 
Diffuse to direct ratio measurements in the 320 to 340 nm region 
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can b« used to Infor •ffoctiv* aerosol optical depths here since 
in this region the ratios are rather insensitive to tropospheric 
OBone. The effective optical depths can be extrapolated to smaller 
wavelengths. These effective optical depths will be somsvdiat dif- 
ferent from the true aerosol optical depths and will tend to com- 
pensate for errors in the ground albedo and aerosol single- 
seatter albedo. The diffuse to direct ratio around 300 to 320 nm 
is insensitive to the altitude distribution of the aerosols so 
long as the ozone profile near the ground is not changing too 
rapidly. Also, detailed knowledge of the stratospheric ozone 
structure or thickness is not required. 

V. CONCLUDING BEMARXS 

Strictly numerical methods of inversion are becoming pre- 
dominant in remote sensing these days. These are, of course, 
valuable to Infer the irregularities and statistical fluctuations 
in atmospheric properties. The analytic model method which we 
have illustrated can be a valuable supplement to such numerical 
methods. They are particuleurly useful when used in conjunction 
with dynamical models of atmosphere structure because of the 
additional physical input of such models. When the models are 
joined to ground-based point sampling data, this remote-sensing- 
analytic model approach gives approximate answers to important 
questions involved in many public policy decisions on atmospheric 
pollution. 


SYMBOL 

a size distribution parameter, corresponds approximately to 

the size of particle where n(r) peaks 
h ozone distribution parameter which is proportional to 

width of ozone density fvinction; prime denotes another 
value of h 

m real part of aerosol refractive index 
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abiolute seattwrad Intaiwlty 

n(r) dlffacantial aerosol size distribution 

M(r) cumulative aerosol size distribution 

II size distribution parameter which is equal to total num- 

o 

bar of aerosol particles 
r aerosol particle radius 

total ozone thickness > prime denotes another value of w^ 
w(y) the ozone thickness function 

y altitude 

y^ altitude at which ozone density peaks; prime denotes 

another value of y. 

o 

ZK zenith angle 

0 scattering angle 

V size distribution pararaetgp; which determines power law 

dependence of n(r) at large r 
p(y) ozone density function 

T optical depth 
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DISCUSSIONS 

D 0 irmnidjimnt With regard to the importanea of multiple ecattering 
in the aureole, I don't want to blow my horn, but Z certainly wish 
to protect Professor Sekera's memory. In my 1956 doctoral thesis, 

I looked at the problem of the solar aureole at his suggestion. 

The method used was essentially Sekera's idea and consisted of 
treating the aureole as a perturbation on the Rayleigh multiple 
scattered skylight. And that is precisely what I did. Have you 
seen my original 1957 paper on the aureole? 

Ormmnt Yes, we have seen your paper. 

Deirmandjian: It was exactly that. You increase the optical 

thickness of the Rayleigh atmosphere by adding a perturbation 
optical thickness due to -the aerosols. When you do that, it is a 
kind of hybrid method where single scattering on the aerosols pro- 
duces the aureole, but multiple scattering, mainly on the Rayleigh 
particles, produces the rest of the background skylight. I thought 
you didn't make that clear. 

* 

Green t We weren't aware of that aspect of ynur paper back then. 

We were aware of your work on the aureole and we have quoted it in 
our work. 

Dairmendjian: 1 think that was the principal point. I have since 

reexamined this method in 1970. I did not publish the results in 
the open literature due to lack of funds. But they are available 
in a formal 1970 Rand Corporation report, in which I introduced 
new phase functions. Indeed, the curves that are obtained look 
very much like some of the measurements I have seen. Subsequently, 
I intended to compare them with the measurements and look into 
their use to get information about the aerosol size distribution. 

At the time, I was unable to do this for lack of support. 

Graani Well, there was no intention to slight you or Professor 
Sekera. We have in our work in this area acknowledged this. The 
UV problem has a new interest ip light of its biological aspects. 

So part of this work was directed toward answering some particular 
questions about the radiance of the UV aureole where some serious 
problems remained. My point about the low altitude ozone distri- 
bution is that you can take advantage of the extra scattering 
associated with aerosols to extend the path and the absorption in 
the low altitude ozone layer. Thus, you can actually drag out a 
little Information about the high altitude distribution, but little 
about the low altitude profile. 

trainman: While I agree with the previous speaker that the aureole 

technique is a powerful one for determining size distributions, I 
should point out that invisible cirrus clouds can plague this 
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tachnique. While there ere frequent casee where you can get nice 
■Mooth functions to which you can apply the theory* this hazard 
always exists. If one looks at this data* one can pick up cirrus 
clouds which are not at all visible to the naked aye. I think 
that it takes a certain amount of judicious discrimination to get 
eases to which you can apply the radiative transfer theory. 

Ormam 1 think Adarsh has shown some results where contours of 
isodensity on a film show distortions that are suggestive of a 
thin cloud layer. In some of our work* we have noted some problems 
of that nature* which one can* perhaps* use as information content. 

Deapakt Yes, we have been looking at the shapes of the solar 
aureole Isophotes rather carefully and we see small bumps. The 
isophote curves arc not very smooth. That could be due to the 
presence of thin clouds. In fact* we have an airport nearby; 
when aircraft take off we can easily detect the presence of con- 
trails from the systematic distortion in the shape of the isophotes 
over the region of the contrail image. 

Fymatf I was very interested in yous conclusion in determining the 
scale height for the troposphere. 

Green t For the tropospheric ozone distribution. 

Fymat: Yes* but you need* as you say, the ozone distribution at 

the lower altitudes. However* it seemed to me, and Dr. Mateer may 
care to comment on it, that the actual ground-based Dobson measure- 
ments are really being phased out in view of the fact that they 
produce different results on inversion using different methods. 

So* if we cannot get the distribution at the lower altitudes, how 
practical is your conclusion? 

Green: Hell, Insofar as profile, the Umkehr method (not the Dobson 
method) relies, and Carl can correct me since he is much more 
expert* on the setting of the sun and inverting sky radiances. Our 
method could be worked at high noon since it is instantaneous. He 
would use wavelength information and we have, in effect* shown by 
simulated tests that the diffuse to direct ratio is only sensitive 
to the addition of extra ozone below about 10 kilometers. So that 
while it is not a complete profile, if you use it in conjunction 
with ground level chemical measurement, it gives you information 
which relates to the lower atmospheric ozone and it is insensitive 
to the stratospheric ozone which is far more abundant. Now am I 
correct that the umkehr method more or less gives you the most 
Information about the higher atmospheric ozone? So there is no 
real contradiction with anything previously known. 
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The unkehr measurement is a different kind of measurement. 
Xt is a ratio measurement of two wavelengths as opposed to what you 
are doing->a ratio of diffuse to direct sunlight. 

Green/ Yes. He also are helped if vm use tw wavelengths as well. 
It eliminates a little bit of the absolute radiometry problem. It 
gives us another check. 
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King I Alaxr what do you think we can laarn about solving the 
atmospheric temperature inversion problem from the similar type 
inversion problem in interring nuclear potential in nuclear phyrics? 

Green; Welli I was discussing. this type of thing with Dr. Wark> 
and the analogy would be if we brought to bear to the problem all 
the physical dynamics data, not only at one location but at adja- 
cent locations. Thus, if we had relatively few questions that we 
posed to the radiative transfer problem, for example, to get hori- 
sontal gradients in temperature profile, and perhaps include such 
facts that then’s tends to be a tropopause and we are mainly inter- 
ested where the temperature break occurs, and ask questions which 
embody all of our experience as to reasonable temperature profiles. 
Then when we apply this external information to an inversion prob- 
lem. it is usually much easier to come to a physically reasonable 
answer. Now l don't know if that's a good paraphrase of the 
analogy. He do have one problem which I have to confess compromises 
the analogy. In the case of the nuclear force problem we think all 
neutrons and protons are the saade, and unfortunately the atmosphere 
takes on so many different states that we dc have that extra com- 
plication. Thus, if you tried to unfold many, many details of the 
instantaneous atmosphere, you may be in trouble with this approach. 
On the other hand, if you are satisfied with answering the types of 
physical questions that, our dynamical meteorologists, pollution 
experts, or biologists want answered. I think you get a lot of use- 
ful information by these modeling approaches. 

Herman: Alex. I would like to ask you. what is the sensitivity of 

the direct diffuse — insofr.r as solving for ozone — to variations in 
aerosol contents? Because as you know we are trying to apply it to 
learn various parameters of the aerosol. Now you are using it to 
learn about the ozone but it still has a sensitivity to aerosols 
.even at the ultraviolet wavelengths. 

Green: Vss. I should, unless I commit another error of omission, 

mention that Ben has used the .direct-to-dif fuse ratio to infer the 
imaginary part of the index of refraction working at 500 nanometers. 
We were not aware of this work when we started ours, but were av/are 
of it during its course. I think the trick in our case is to 
choose two wavelengths — one in the region where the ozone is 
absorbing ani one where it is nonabsorbing. Then we take advantage 
of whatever we know about our aerosols, including a ground-level 
aerosol measuring device, plus an aureole measuring device. I 
would not go to anything expensive like a lidar. But I think if 
we bound our aerosol model somewhat we find by sensitivity analyses 
that we do get some information about the low level ozone. This is 
what we would pick up by using two wavelengths, just as Dobson does. 
We would go to about 320 nanometers and then to 305 nanometers. We 
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would oxpeet that tho ••rosoln do not chango vary Mieh« ao that wo 
hold tho aoroaola eonfutlon factor do«m. 

GouJardt Z an trying to apply Invoraion tochnifpioa to tho flold 
of conbuotion. Z got tho piofiloo of tonporaturo or concontrationo 
following tho nethodo that Dr. Chahino hao dovolopod. But. Z got 
than in toma of optical thicknoaa. And whon Z want to convort it 
to phyaical thicknoaa. which io what tho ccmbuation peoplo want. Z 
find that Z cannot depend on the conatant concentration of CO 2 aa 
you do in the atamaphere. Could someone give me acme help on how 
to get around thia problem in frequency acanning only? 

Chahino f Iiet mo aitawer Profoaaor Goulard 'a quoation on tho aolu- 
tion in torma of tho optical dopth. You know, the independent 
variable in the radiative tranafer equation ia t but we aet dx ■ 
Ox/lxlda in the equation and obtain a aolution aa a function of 
tho phyaical acale z. Aa you know dx/3z ia obtained on the baaia of 
an atmoapheric model. Thua. the aolution ia truly a function of 
T although it ia uaually preaented aa a function of z. The trana- 
formation f roai the t-acale to the ^z-acale ia done through a model . 

Goulard t You aaaume Laplace Law of the atmoaphere and a conatant 
mixing ratio? 

Chahinet Yea. we have to do this. 

Goulardt And for a flame? 

Chahinet You would have to develop a model for your flame and 
correlate your optical depths with a physical scale. 

Goulardt la the concentration of CO 2 really constant throughout 
the atmosphere? 

Graent Hell, it ia. yes. within a small period of time. There 
are some people who think it is growing. 

Cilia t The concentration of CO 2 in the troposphere is ^ite con- 
atant although not absolutely so. It varies seasonally by perhaps 
one part per million out of a background of 330 with a long-term 
trend of. perhaps. 0.7 ppm per year. 
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^ OF ultraviolet SKYLIGHT 

A. Ej S. Crken, K. R. Cr(>» and L A. Smitii ... 
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AkMrar*— We present nn improved analytic characterization of diffuse spectral irradiance (slcylii;ht| for 
the wavekn>;tH range SSO-.'tSOnm and sislar /enith angle raii;v from 0 to 85'. The formulas achieve 
treater accuracy by (a| fiKusing on ratio represenlatio'AS and (h) adjusting the purnmeiers to the more 
precise radtaiive transfer calculalions of Dave Dra^lau and Halpern. 


INTRODUCriON 

There is a serious need for simple but accurate analy- 
tic characterizations of skylii>hl, particularly in the 
' UV where skylight soinetime.s has greater biological 
importance than direct sunlight. Previous work by 
Green tt at. (1974) (GSS), Shcttic and Green (1974), 
Green and Mo (1974) (GM), Mo et al. (I97S) and 
Green (1976) have provided simple approximate 
phenomenological representation of the systemutics of 
Bener (1972) and of numerically calculated diffuse 
. spectral irradianccs. However, experimental studies of 
Chat and Green (1976), and Garrison et al. (1978a, b) 
have focused atti vtion on the ditfuse/direct ratio 
which avoids problems with absolute radiometry 
which are particularly difRcult in the UV region. 
Furthermore, this ratio depends more slowly upon 
wavelength than the diffuse spectral irradiance itself. 
Unfortunately the dependence upon angle of the dif- 
. fusc/direct ratio can be large because the denominator 
(the direct irradiance) becomes so small at large solar 
senith angles. In the present work we will use more 
ehborate ratio techniques to find an accurate analytic 
characterization of skylight which^ for all angles and 
wavelengths of interest, will bring the At to numeri- 
cally generated diffuse irradiance from the 2S% level 
of accuracy to about the 5% level. Toward this aim, 
rather than use the numerical output of the code 
available to us (Slicttle and Green 1974) we will use 
the results of the more precise radiative transfer calcu- 
lations of Braslau and Dave (1973) and Dave and 
Hatpern (1975). Their results (to be referred to as 
DBH) for six models (B, C, O, Cl, Dl and Cl*, sec 
Table l)have kindly been made available to us in the 
fortn of magnetic tape output and readable hard copy. 
In clfcct, the objective of this paper is to make the 
results for the various DBH models available in ana- 
lytic form fur biological users. The form developed is 
more elaborate than those previously proposed but 
should be worlluvltile if n precise knowledge of tliu 
State of the atmosphere is available. 

ANAt.VriC FOUNt FOU IMRKtT 
IKMAOUNCK 


The problem we arc involvctl with has many' 


dimensions, X, the vv, -rgih, 0. the sun angle, wj, 
the ozone thickness, Wj. he aerosol thickness, (both 
of which arc dependent upon y, the altitude) and A, 
the ground albedo. In addition, the proAlc dwracler- 
istics of the ozone and aerosol enter. Rather than 
begin at the beginning in this attempt to describe UV 
skylight as a function of these many degrees of free- 
dom, we shall attempt to utilize the experience gained 
in our earlier studies. 

The downward* direct spectral irradiance at the 
ground may be placed in the basic form. ] 

D{X, 0) » cos 9 //(;.)cx‘p - 1 (t,/#ij) (I) ! 

1 . I 

where T|, Tj and tj denotes the air, aerosol and ozone | 
optical depths, respectively, and /<, jq, jij and /tjj 
denote cosfl and generalized cosine functions which 
arc appropriate to the three species in view of tire 
roundness of the earth. We will express these func- 
tions in the form 





( 2 ) 


where the (,’s are smalt characteristic numbers which 
depend upon the altitude distribution of the species. 
(See GSS 1974 for a more complete discussion.) 


atARACTERIZATION OF SKYLICIIT 

It would be time and space consuming to recount ! 
all the analytic forms, mathematical models and | 
approaches we have examined during the course of 
this investigation. In addition to molecular band type I 


Table 1. Atmospheric niudels 
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Gaseous 

Height 

Refractive 

Model 

absorption 

distribution 

index 

A 

No 

No 


B 

Yes 

No 
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C 

Yes 

Average 

1.5-O.Oi 

D 

Yes 

Heavy 

t.5 0.0i 

CT 

Yes 

Average 

1.5 OOli 

Dl 

Yes 

Heavy 

I.S-OOli 

Cl* 

Yes 

Average 

1.5 0.011 
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Figure I. Showi our fits to i^(A,0). The smooth curves are 
■ our analytic fits and the points reprsiicnt the ‘data*. 


models we tried again to adapt the two stream model 
(Mo ft o{., 1975) for the purpose at hand. However, 
'inaccuracy or cuihbcrsomcncss led us to reject many 
vpproaches in favor of the analytical model which 
finally evolved. In this model our basic strategy is to 
represent the downward dilTuse spectral irradiance (S) 
at the ground in terms of two ratios and the direct 
spectral irradiance for an overhead sun (0 0 , or 

|i >■ I). Thus fpr dilTusc sky irradiance we use 


5(40 - y(;.,0)A/(A)H(A)cxp 


-(?") 


(3) 


where 

with 


y(A.fl) - S(A.(n/S(;..0*) 


(4) 1 

(5) 


Thus, representniions at the tew per cent level of 
accuracy should be feasible except possibly at the lar* 
gest solar zeiiilli angles and the shurlest wavelengths. 

Tlie points in Pig. I illustrate variatkitis of .V with 
0 at H'j m 0l 30 using 'data* from Dave and llalpern 
(1975) Model C. Here the points are numerically | 
generated results at 297.5. ^>7.5 and 330.0 nm. Tiiel 
smooth curves are derived from an analytic rcf resen* ! 
tatkm in the form ) 

0) • [F + (I - F)cxp - yj(t, + f 4 )tMexp I 

■^■(yi»i + yi*i)^ («) 


where 


1 


I + t4iV»(y) + 0fcs(*)r’ 




g-go + a,y, (7) 

+ (*) 


whem here T 4 denotes the part'iculate absorption opti- 
cal depth corresponding to Tj the particulate scatter- 
ing optical depth, A-j(A) the .ozone absorption cocIFi- 
cicnt.'Af|(y). Afj(y); and Siiy) are the total thicknesses 
above the altitude iv) each normalized t.o unity at sea 
level (y B 0). The mathematical form we use for r.'ch 
of these is (see Fig. 2'and Table 2). . 


X 


N.(y)-/t 


1 + «i 


Bi + exp0’//ii) 


The virtue of this combination of ratios is that the ‘ 
dynamic range of M is less than one order of nuigni - 1 
tude for the wavelength range of intcrc.st (2itO-3HO nm) 
and for various aerosol optical depths and ozone ^ 
thicknesses of interest. Furthermoie, the dynamic | 
range of is also reasonable (I onlcr of magnitude; 
above 320 nm, 2 above JtXtmn. 3 orders above* 
290 mn nnd 4 orders of magnitude above 2S0nm). j 


..... »+exp-yr/X. 

where yi denotes the altitude of the maximum density 
of species i, /it and yt are scale height parameters, Xf is 
a shape parameter, and/j is a fitting parameter. All of 
the parameters in Eq. 9 arc dependent on the model 
distribution used. The first tcrmJi^Eq. 9 controls the 
low altitude distribution. thVsccond'the 'high altitude 
distribution. With this representation the T|(A) is the 
optical depth at sea level and ti(A)/y/| 0 ') is the optical 
d^th above any altitude. Tlic Rayleigh and ozone 
extinction coefficients used by DItH are given in 
Table 3. Tlic spectral sha|ie of the aerosol extinction 
coefficient over the range 280-350 nm depends upon 
the size distribution assumed. For the DBH models it 
appears adequate to use an average aerosol optical 
depth independent of wavelength. Table 4 gives the 
best fit parameters obtained for tlic UIIH models. 

Next wc examine the ratio Af. After trying a 
numlier of approaches, we obtained favorable results 
for the DBII models with the representation 

Af(A) 


-4.it, W,(.v) + /l.j[r,A',Cv)r 


1 + /•.jHjiV,(y)[r,N,tv) + uNtiy)}"- 


( 10 ) 


Table 4 gives the Af paiamctcrs for the DBII . 
models. Figure .ta illustraies the fits as a function of | 
wavclengtli for Model Cl* with wj « 0.2, 0..3 and 0.4. j 
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Table 1 Adjusted parameters for Nfyl 


V 

• . . t 

• . . I • 

‘ Average aerosols' 
2(C) 

Heavy aerosols 
3 2(D) 

/, 

1.000 

0.960 

a039 

0.956 

dt 

0A31 

-0.I4S 

2350 

-0.448 

•*. 

6.350 

0.952 

. 2.660 

1.003 

Ft 

— 

16.330 

22.510 

19.378 

Zt 


3.090 

4.920 

1.106 


1.8 X I0-* 

3.0 X 10** 

lA X I0-* 

33) X 10*‘ 


Table 3. Wavelength dependent optical parameters as used in computations 


1 

Wu'AslengtIi 

(pm) 

m 

Solar energy 
incident at tup 
(mWcni‘*/<m' *) 

Nominal 

s|>cctral 

bandwidth 

l/im) 

t,(2) 

Scattering 
optical 
thickness 
due to 
molecules 

m) 

Ozone 
absorption 
coellicicnt 
(atni-cm)* ' 

T,(A) 

Scattering 
optical 
thickness 
due to 
aerosols 

t«(7.) 

Absorption 
optical 
thickness 
due to 
aerosols 

0 : 297 s 

54.90 

0.005 

1.2652 

13.58 

0.08012 

0.01369 

03025 

5S.SS 

0.005 

1.1779 

6.76 

0.080.52 

0.01355 

03075 

64.60 

0.005 

1.0980 

3.50 

0.08091 

0.01.342 

03125 

72.65 

0.005 

1.0248 

1.73 

0.08129 

0.01328 

03175 

79.70 

0.005 

0.9.577 

0.890 

0.08166 

0.01.315 

0.3225 

90.35 

0.005 

0.8960 

0.460 

0.08203 

0.01303 

03300 

105.90 

0.010 

0.8123 

0.1730 

0.08256 

0.01282 

03400 

107.40 

aolo 

0.7164 

0.0515 

0.08324 

0.01257 

03500 

109.30 

0.010 

0.6341 

0.00S6 

0.08.390 

0.01233 

0.3600 

106.80 

0.010 

0.5634 

0.0014 

0.08452 

0.01209 

03750 

115.70 

ao20 

0.4751 

0.0000 

0.08539 
f, -aiKSI* 

0.01175 
I 4 - 0.01 3* 


* For Dl Moilcl multiply t| uml (4 by 4.18. 
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Table 4a. Paramcleri for y, M ami r(i, h j) for average mode! 




Species 

Rayleigh 

Particulate 

Ozone 

parameter 

Im 1 

1-2 

1-3 • 

A 

0l5777 • 

a427 

ft977 

Ail 

0.7H79 

I95S0 

ai97t 

•• • Ati . 

ftS399, 

0i023f> 

E67I0. 

: A„ 


* 

3.2SS 

a 

Table 4b. 


Alpha 

parameter 

SbQ 

a»a 

a ■ b 

<1. 

1.104 

1.0790 

0.SV66 

Pi 

— 

I.S330 

0.7847 

fi 

1.433 

0.067 


ta 

aoi99 

-3.31 X 10** 



GROUND ALUEOO DEPENDENCE 

It is possible to infer the global irradiance G(A, i, 0) 
for any ground albedo (/I) from the global irradiance 
C(Q> A.0) for zero ground albedo by using the formula 


C^A.0)> 


r- r(;.M 


01) 


where r(A) are the air reflectivities. The air reflectivity ; 
values of Shcttic and Grech have fit rather accurately i 
by a similar equation ns the analytic equation for : 
Af(A) above i.e. 


from the air and aerosols would not strictly by Lain- 
bertian. However, using to correct the downward 
diffuse spectral irradiance to the equivalent ground 
level but for zero ground albedo we use 

S(0, y,) - S(0: y) - » y)U(0./. y) . (14) 

where r().,'y) is the air reflcctivily fonction. We would » 
expect this function to fall ofl' with altitude in a wuy 
dependent upon the air and aerosols . above the 1 
ground. Since the backscaltcr of air is usually much ! 
greater than the backscattcr of aerosols, we. will let i 




+ /•wCTsWify)]'’* 


( 12 ): 


I + A„[T,^,(y) + T^V,0)r 

Table 4 gives the corresponding parameters for the 
various DDH models. Figure .^b illustrates the fits for 
the models at wj » 0.2, 0.3 and 0.4. 

The dependence of the diffuse spectral irradiance 
upon the height of the ground above sea level has 
been a source of confusion. Strictly speaking, one can- 
not use the downward diffuse spectral irradiance at an 
altitude y above the ground (at sea level) as equivalent 
to the dow nward spectral irradiance at ground which 
b at a height y, above sea level. Yet this is tlm oiiiy' 
data available to us at this time. We can, however, usc| 
such data to fill out our analytic representation in the 
following way. j 

First, we calculate the equivalent ground albedo, 
associated with the air and aerosols below ground' 
altitude y by using the altitude dependent irradiance | 
S. D and If in I 


’ w 

ot 

ftS 

or 

os 

os 

04 


.02 


DM + S(y) 


( 13 ) 


where V is the upward dilfuse irradiance when the sea 
level ground alltedo vanishes. The idhedos so 
obtained me slightly dependent upon the sun angle as 
' might Ik! expected from the fact that the bad scatter j 
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l•■iBure 3a. Shows our fils to A/|f.). Ag.iin, the smooth 
curves are the fils and liic |it>iiUs are the data. 

I'igucc 31). Shows anal) lie fils lu r(3., 1 V 3 ). 
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height above scsi level according to 

the approximation 

r(;,y) « r(^..0)Ky)/Po “ r(;, 0 ).V, 0 -). (15) 

This permits us to generate diffuse spectral irradiance 
values for ground level at an cquivuteni zero albedo. 
We next use the procedures iu the previous section to 
analytically represent We have tested this 

procedure with model C. The altitude dependent coef- 
ficients are well characterized by Eq. 12. 


WAVtXENCTII DEPENDKNT FUNCTIONS 

Previously Green et at. characterized the extrater- 
restrial solar spectral irradiance ll{y) by 


Tlie Rayleigh optical depth can also be fit reason- 
ably well ill this s|icclial region by 


m) 




(16) 


This provides a reasonable representation in the inter- ' 
val 290 to 3‘tOnm. However, to extend this we will use 

«(;.) - /cjl - cxp[ - K- exp^i^j jj (1 7) j 

where K » |^>9.5. k = 0.6902 , 6 => 2.7.74, - 300 nin. ! 

Figure 4 shows the fils in relationship to the data of 
Thckackaia (1971). It should be noted that the Planck ' 
function docs not lit this spectral region no mailer 
WhitI temperature is chosen, .. ! 


where 


»sW-t,o(^Y*’ 


Tif ■» '1.221, 


. <I8) 

' Figure 4 also shows this behavior. 

{ Finally the ozone attenuation cocITicient may be 
represented quite well by 




^ + L 


I 


fi + expl 


C-^) 


(19). 


7.294. 


where ko "^•5I7» fi “ 0.0445 and 6 

Figure 4 also shows this function. I 

For the aerosol models chosen by Dave et at. it is | 
sufficient. to treat Tj ns a constant, i.c. its value at' 
b 300 nm. For other models we could revert to the . 
expression used earlier by CSS, i.e.: ; 

t,-t,o(VA)** (20)1 


. OISCySSION AND CO.NCLUSION 

Wo have developed a mathematical expression for 
the diffuse spectral irradiance (skylight!) which has 
been adjusted simultaneously to thp various models of 
DBH using their extinction and attenuation coefii- 
cicnls. The fits to these data arc generally within \V,i 
for the wavelengths (A) zenith angles (0), ozone thick- 
ness (wj) and aerosol optical depths (r.) available as 
input data. The. fits have been generalized to encom- 
pass several ground albedo (A) and ground altitude 
(y,). The function of six variable (A, 0, W 3 , tj. A, yjt so 
obtained can now reasonably be used for various in- 
terpolations and extrapolations between or somewhat 
beyond the data points used in adjustment process. 

Our equations are capable of more precise fits to 
numerical data (5%) if one limits the number of 
models that are fit simultaneously and if one uses 
numerical inputs for //(A), k)(A) etc. Thus we can de- 
scribe the Cl* model output, for example, with con- 
siderably greater precision ami wun fewer parameters 
than we <xtn by fitting all the models simultaneously. 
Table 5 contains the parameter values for the indivi- 
dual fits to the various models. For some specialized 
purposes such representations might be useful as 
alternatives to the original table.s. 

For the purposes of most of the applications our 
generalized nuulel adjusted to all of tlic cloudless 
DUIl models should be most useful, l-'urlhcrmorc, if 
we replace the wavelengih dcpimdent in|)uls used by 
DIUI by the smooth mathematical functions de- 
scribed in the Wavelength Dependent Functions sec- 
tion we can coiiipute approxintate spectral irradiances 
at arbitrary wavelength separation. 'I'his shouhl be 
lielpful to photitbiologists who ntust convolute spec- 
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Tables. 


Model 

Parameter 

D 

C 

Cl 

€!♦ 

D 

Dl 

(a) Af(A) Parameters 




f 


^•1 

0.7120 

0.M69 

0.7994 

0.862K 

0.8012 

0.7994 

^•1 

— 

-0.7872ie«)<.U274 14.614 

3.6567 

3.3274 

^•1 

0i2»} 

03016 

UI9IK 

0.2S25 

0.1556 

01911 

P. 


3.2873 

1.6552 

7.7010 

1.6368 

3.2873 


l.0.^S8 

0.8909 

I.0C9I 

0.9257 

1.0674 

1.0691 

(b) rii, w«) Parameters 






^¥1 

O.S332 

0.5640 

0.5595 

0.5608 

0.4789 

0.4375 



OlOOO 

-0.0I27»«}<0.I068 

0.1470 

0.07316 


0.SK08 

0.6470 

0.6543 

0.6428 

0.5197 

0.5829 



1.0000 

0.81 i8 

0.2455 

0.5165 

0.34.30 

«k 

0.6S39 

0.6040 

0.6020 

0.6050 

0.6334 

0.6102 

(c) £^U, 9) Parameters 






Yi 

0.542 

0.5842 

0.5565 

0.4758 

0.5338 

0.4840 

It 

0.U 

0.3968 

0.5583 

1.1383 

0.6832 

0.7779 

U 

a9833 

0.9968 

1.0046 

1.0158 

1.0234 

1.0184 

o 

1.0712 

1.1826 

1.2010 

1.0210 

I.ISIO 

1.2083 


I.0S84 

7.7196 

I.0II7 

19.98 

4.7746 

0.63SS 

f. 

0.01992 

0.02285 

0.02091 

0.02253 

0.02892 

0.03825 

9 , 

0.S94I 

0.4123 

3.0985 

a270l 

0.5789 

3.4702 

9 * 

0 

0 . 

0 

0 

0 

0 

t. 

0 

0 

0 

0 

0 

0 


I 


< 


I 

i 

I 

I 
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Table 6. Analytic models of ^(A, 0) tested* | 


Model ' 

* ’ Formula 

Fit 

Exponential 

e‘*4 

Poor 

Double Exponential 

(!•/>-•♦ +'Fe-»* 

. Best 

Mayer Goody ^ 

f >iMI * Mt 

Good ■' 

Verhulst ^ 

FHe>* - 1 + F) 

Fair 

Vcrhulsl 0 

(p + I)/(P + e*») 

Poor 

Verhulst 0* 

(P + l)/iP + c'*’) 

Fair 

Mayer Goody 0^ 

Q-itVI ♦ l»»lt 

Fair 

Modified M-C 0^ 


Fair 

Gaussian -1 

e”*** 

Poor 

Gaussian -2 


Good 

Gaussian -3 

(1- F)*-’»’ + Fc-4** 

Fair 

Fraction to Power 

[a/(l - /I + «)]* 

Fair 

Malkmus 

exp - 2a//lt(l + - 1] 

Good 


' * ^ is denned by Eq. 8. The symbols fl,y,S, and f arc wavelength - 
dependent parameters. 


Ir 


tral irradiances witit various types of action . ‘tcctra to ; 
arrive at an equivalent dose rate for studies of dose | 
response relationships. For such purposes the rrcgii- ! 
lar structures in //(/), k^iA) should uverugu out dr the , 
most part. j 

Our analytic altitude distributions can also serve as I 
a convenient way of inputting other possible atnios* ' 
pneric models. It is not unreasonable to expect that ' 
the spectral irradiances for somewhat varied altitude | 
distributions will change in a fashion called for by our [ 
equations. It is also reasonable to allow for other ; 
absorbing components such as hii'.h SO^ concen* ! 
Irutions in aq ui km area into our model by incorpor* ; 
ating an additional w^ki(A)i\i{y) in parallel with the 
ozone term. I 

It might serve a useful purpose to inventorize thel 
alternatives which we have examined during the; 


course of this study. The GSS formula and the GM| 
variation were first considered but set aside in favor! 
of models which led to simpler expressions for thej 
ratio (Fq. <1.) The two stream formula by Me 

el at. (1974) was also considered but was found to be! 
very cumbersome when modified to achieve uccuracy.| 
Table 6 lists .some of the formulas which we have; 
considered for the y{0,A) along with comments I 
as to our fitting experience. While wc cannot say that | 
further data fitting will not change the overall rating | 
the double e.xpoiienlial has the advantage of lloxibility I 
and simplicity. The ratios Af(7.) (I'.(|. lU) and r(A, Wj) 
(Fq. 12) were directly adapted from the CSS elTort. 

In final conclusion we trust that in addition to its 
liscruliiess to phntohiologisls our equal ion will serve 
as a frame of reference for charting experimental 
observations ofclilfusu tIV .spectral irradiances. 


I 




CkwactwiMlion o( UV ikylight jf 

Aetmwtntffetmutii—The writcrt wuuld like to thank Dr. ‘ Tliii work was supiKtrled in part by the Naliotial Aer* I 
Norman Hraslau ami l>r. J. *' Dave for making available , onaulics ami Space AdininisiraliiNi Uaperinieni Dclinilion I 


tapes ami paper copy output of their original radiative and Special I'rojccis llianch (under Contract No. I 
Iransfer calculatiuns. We would also like to thank Drs. ' NAS9'I5II4) and the Nimbus 7 Data Application ODicc 
Martyn Caldwell. Kobert Rundel and Richard Zepp for \ (under Contract No. ATM 7S'2l962k 
their helpful comments. j , 


, RETERENCES. .. ... 

Bcncr, P. (1972) European Research Office. U.S. Army, London, Contract No. DAJA37>fi8*C*IOI7, p. 
S9. 

Braslau, N. and I. V. Dave (1973) Clfeet of aerosols on the transfer of solar energy through realistic 
model atmospheres. Part III; Ground level Iluxes in the birdogically active hands, l3Si(y-.37UO mic* 
rons, lUM Research Report. RC 

Dave. j. V. and P. Il.-ilpcrn (1976) /Itmor. Cna. 10, S47-SSS. 

Chai, A. T. and A. E. S. Green (1976) Appl. Opt. IS. 1 183-1 187. 

Garrison. L. M. L. E. Murray. D. D. Doda and A. U. S. Green (1978a) Appt. Opi. 17, 827-836. 

Garrison, L. M.. L. E. Murray and A. E. S. Green (l97Sb) Appl. Opt. 17, 683. 

Green, A. E. S. (1976) Solar spectral irradiance reaching the ground. In ProMdings on Nonbiolofiical 
Tnmport and Ttam^ormation of Pollutants on Loud and in iVater. National Bureau of Standards, New 
York. 

Green, A. E. S., T. Sawada and E. P. Shettle (1974) Photochtm. PhotohioL 19, 33 1 -359. 

Green, A. E. S. and T. Mo (1974) Proceedings of the 3rd Conference on ClAP, DOT-TSC- 
OST-74-15, 518-522. 

Halpern P., J. V. Dave and N. Braslau (1974) , Science 186, 1207. 

Mo, Th G. P. Kerwer and A. E. S. Green (1975) J. Ceophvs. Res. 80. 2672-3676. 

Shettle, E. P. and A. E. S. Green (1974) /Ippf. Opr. 13, 1567- 1 SSI. 

Thckaekara, M. P. (1971) Solar Electronuignetic Radiation. NASA SP-800S. 





AEROSOL EFFECTS ON ATMOSPHERIC RADIATION 


( . 7 ‘? 

6 


IN THE MIDDLE ULTRAVIOLET 

A. E. S. Green and P. F. Schlppnick 
University of Floridafl 

The needs for an accurate characterisation of ultraviolet (UV) skylight 
in several disciplines are sometimes better served by approximate analytic 
representations rather than by large masses of computer outputs. Toaards 
this end we first review sources of numerical "data” for UV atmospheric 
radiation and prior attempts to characterise UV skylight analytically. 

After indicating the various formulas examined we present illustrations of 
fits of the Green^ Cross and Smith formula to theoretical spectral irradi- 
ates calculated by Braslau and Dave. The degrees of fref>dom include wave- 
length, solar senith angle, osone thickness, aerosol content, ground albedo 
and altitude. Finally, we present a representation of Peterson's actinic 
UV flux in a form similar to that developed for spectral irradiance. 


To be published In the Proceedings of 
Atmospheric Aerosols Workshop 
November 6-8, 1979 

^ICAAS, Gainesville, Florida S2611 


1. Introduction 


Accurate characterizations of ultraviolet (UV) skylight are needed 
in several disciplines. These Include photobiology* where UV plays an 
important role In vitamin D production* In skin cancer dose response 
studies and in agriculture where It adversely influences plant growth* and 
in biological oceanography where It Influences phytoplankton. In addition* 
UV impacts upon air pollution* tropospheric chemistry in general and in- 
directly upon climate change questions. The need for accurate characteri- 
zations Is probably better served by approximate analytic representations 
rather than by large masses of computer outputs. Towards this end we first 
review sources of numerical "data" on UV atmospheric radiation. Then we 
review a number of attempts by the Florida group to characterize UV skylight 
analytically. These began In 1973 with the Green* Sawada* and Shettle 
representation and continue with the 1979 formulation of Green* Cross and 
Smith. After Indicating the various formulas examined we present Illus- 
trations of tits of these formulas to experimental and theoretical spectral 
Irradlances. The degrees of freedom Include wavelength* solar zenith 
angle, ozone thickness* aerosol content* ground albedo and altitude. 

Finally* we present a representation of actinic UV flux In a form similar 
to that developed for spectral Irradlance. The purpose Is to provide 
the user with convenient access to the massive body of radiative transfer 
calculations in the middle ultraviolet regions. 

2. Prior Work and Sources of Data 

Experimental and theoretical knowledge of atmospheric spectral 
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irradlanc* In the ultraviolet (UV)* particularly the dlffuae component 
(akyllght) waa at a very primitive atate until the 1960a. From 1958 
to 1969 Bener (1963* 1970a* b) conducted extensive series of measure- 
isents of solar ultraviolet (UV) radiation In Switzerland* idilch to this 
day provides the best available body of spectral data. The systematica 
of these measurements as functions of solar angle* wavelength* ozone 
thickness and other variables were gathered together In a final report 
(Bener* 1972). 

The first major theoretical attempt to present UV spectral Irra- 
dlances reaching the ground was by Dave and Furakawa (1966). They pre- 
sented solutions to the radiative transfer equation for Rayleigh atmos- 
phere having realistic levels of ozone. Serious attempts to allow for 
aerosol scattering came out first In 1973 In several parallel and Inde- 
pendent efforts. Green* Sawada and Shettle (GSS* 1973* 1974) used a 
phenomenological modification of a formula based upon single Rayleigh 
and aerosol scattering adjusted to Bener *s systematics as an analytic 
representation of the diffuse spectral Irradlance (skylight) . Simul- 
taneously Shettle and Green (1974) developed a fast multiple scattering 
calculatlonal method to predict UV spectral Irradlance. Shotkln and 
Thompson also In 1973 used atmospheric models containing aerosols 
along with a Monte Carlo method and the two stream approximation to arrive 
at UV spectral Irradlances. They compared their results with Bener *s 
data and found some discrepancies of the order of 20%. Braslau and Dave 
(1973) also about the same time* used direct solutions of .the spherical 
harmonic approximation to the basic radiative transfer equation to 
arrive at diffuse UV spectral irradlances. Thus by 1975 our knowledge 
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of UV radiation levels and Che Influence of aerosols upon it had 
advanced very greatly. 

In connection with the applications mentioned in the introduction, 
there is a serious need for simple but reliable analytic characteri- 
zations of ultraviolet skylight which usually have greater biological 

importance than direct UV sunlight. The work along these lines by GSS 

Miller (1974) 

and SG, by Green, Mo, and yy • Mo, Kezwer and Green (1975), and Green 
(1976) provide analytic representations of the systematlcs of Bener's 
1972 work and of numerically calculated diffuse spectral irradlances of 
Shettle and Green. To overcome problems of assigning absolute Irradlances, 
Chai and Green (1976), and Garrison et al . (1978a, b) focused attention 
on diffuse/direct ratios. These ratios depend more slowly upon wave- 
length than does the diffuse spectral Irradlance itself. Unfortunately, 
because the direct Irradlance becomes so small at large solar zenith 
angles, the dependence upon angle of the diffuse/direct ratio is very 
large. Green, Cross and Smith (1979) have recently found a more accurate 
analytic characterization of skylight which, for all angles and wave- 
lengths of Interest, will bring the fit to numerically generated diffuse 
Irradlance from the 25% level of accuracy to about the 5% level. As 
input data they used the results of the more precise radiative transfer 
calculations of Braslau and Dave (1973) and Dave and Halpern (1975). 

This work summarizes the GCS effort particularly as it related to the 
Influence of aerosols on skylight. In addition it describes some initial 
results of an attempt to represent actinic flux with a formula of the 
GCS type. 


r 
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3. The CCS RepresentatT.on of Diffuse Spectral Irradlance 

The spectral Irradlance problem has many dimensions; X, the wave- 
length* e» the sun angle, w^, the ozone thickness, aerosol 

thickness (both of which are dependent upon y, the altitude) and A, the 
ground albedo. In addition, the profile characteristics of the ozone 
and aerosol play a role. 

The direct downward spectral Irradlance reading the ground follows in 
a straightforward way using Beer's Law. It may be written in the form 


D(X,6) « COS0 H(X)exp - 



( 1 ) 


where t^, T 2 * and denotes the Rayleigh scattering, aerosol scattering, 
ozone absorbing and aerosol absorbing optical depths, and (i, ^2 

denote cos6 and generalized cosine functions of the form 

2 ^ 

Mi - (y>t^)/(l+t^) (2) 


Here the t^'s are small characteristic numbers depending upon the altitude 
distribution of the species which allow approximately for the roundness 
of the earth. It would be time- and space-consuming to recount all the 
analytic forms, mathematical models and approaches which have been tried 
to analytically characterize the diffuse spectral Irradlance. The most 
natural model to start with Is the two stream model and attempts to 
use it have been made by Shotkln and Thompson (1973) and by Mo, Kezwer 
and Green (1975). However, inaccuracies or cumbersomeness have thus 
far favored empirical models over the two stream model. In the GCS 
model the basic strategy was to represent the downward diffuse spectral 
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irradiancc (S) at the ground In terms of two ratios and the direct 
spectral Irradiance for an overhead sun (6 ■ 0« or M " 1). Thus 
diffuse aky Irradiance is expressed In the form 

S(A,6) -^a,0)/l{X)II(X)cxp (3) 

where 

J (X,0) » S(X,0»)/S(X,O«) (d) 

with 


;K(X) - S(X,0»)/D(X,0») (5) 

The vlrt'.'> of this combination of ratios is that the dynamic range of 

i^(X,9) and particularly over the wavelength range of Interest 

(280-380nm) and for various aerosol optical depths and ozone thicknesses 

of Interest are relatively modest so that analytic representations at 

the few percent level of accuracy should be feasible except possibly at 

the largest solar zenith angles and the shortest wavelengths. 

The numerically generated points in Figure 1 illustrate variations 

of *^with 0 at w ■ 0.30 using "data" obtained from Dave and Halpem 
oz 

(1975). The dark smooch curves on Figure 1 are derived from the GCS 
analytic representation 

J(X,6) = [ F + (1- F )exp~Y3(T3+T4)4iJexp-(YjL‘r2^ + (6) 

where 

♦ - 1(1 + t)/(M^ + t)]** - 1 


(7) 
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Here denotee the particulate abaorptlon optical d«*pth correapondlng 
to T 2 the particulate acatterlng optical depth and 

P - ^ (I 

1 + AQ3lt3N2(y)+T^N2(y)+g kjCMl^o 

with 


t - t + t y (9) 

o a' 

and 

g ■ gp + g^y (10) 


The functions Nj^(y), NjCy)! and N^^y) are the total thicknesses above 
the altitude each normalized to unity at sea level (y ~ 0). Following 
GSSi GCS and other works they are represented by the mathematical form 


N^Cy) " o^+exp (y/h^) 




l+exp - y^/Xi 
1+exp (y-y^)/Xj^ 


( 11 ) 


With these representations the t^(X) are the optical depths at sea level. 
Also shown on Figure 1 by light lines are the direct ratios 

20 (XfO) - D(X.0)/D(X,O“) (12) 

Note that these fall off more rapidly at large solar zenith angles than 
do the corresponding .^(X, 6) . Thus* as the sun sets the diffuse Irradlance 
Increasingly becomes more Important as compared to the direct Irradlance. 

For 6 ■ 0 the relative magnitudes of the direct and diffuse 

spectral components are fixed by theA|(X) function. The lower curve in 
Fig. 2 gives the/f(» function for w^ ■ 0.3cm. 

l+A^3W3N3(y) [13113 (y) (y) J^a 


( 13 ) 
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Tabic 1 gives parameters for the GCS diffuse spectral Irradiance 
formula for an average atmospheric model. The quantities in parentheses 
will be described In the next section. To complete an analytic repre- 
sentation It Is necessary to specify H(X)» t 2 (X)» and 

T^(X). The wavelength dependence of H(X) In the Interval 280-380mn 
may be represented as a smooth function 


€ X - X N 

H(X) - K ^1 - exp[- K exp j 


<14) 


with K - 3J395 V/ttTnm, < ■ 0.6902, L ■ 23.74, X^ - 300nm. This empirical 
function fits much better than the Planck function no matter what tempera- 
ture is chosen. 

The Rayleigh optical depth and the ozone absorption coefficient 
may be represented by 


T,(X) - 1.^21(X /X)^*^^ 
X o 


(15) 


and 

X - X 

kjCX) » T 2 (X)/w 2 « k^(8+l) / [8+exp( — g — ^)] (16) 

with k - 9.517 (atm. cm)~^, 8 - 0.0445 and 6 - 7.294. For most purposes 
o 

it Is sufficient to characterize the aerosol scattering depth T 2 and 
the aerosol absorption depth as Independent of wavelength over the 
range 280nm to 380nm. While some variation would be expected, the 
usually uncertain differences between atmospheric models from the 
no-aerosol case T 2 * * 0 to an average aerosol case, T 2 ** 0.08, 

■ 0.012, to a heavy aerosol case, T 2 0.34, w 0.05, far outweigh 
these minor wavelength variations. 
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nie foregulng equations apply to the case of zero ground albedo. 

One can use the methodology of Shettle and Green* and Green* Cross and 
Smith to correct for any ground albedo. This entails Introducing an air 
reflectivity function* r(X*v^) which may be represented by 

(17) 


r(X*Wj) - 


^ih^i(y) ^ 

1 + Aj^3lT3N3(y) + x^N2(y)l% 


The complete model so constructed using the average parameters given In 
Table 1 provides a reasonably accurate representation of spectral Irradlance 
for all of the Dave Braslau atmospheric models. 

To correct the Irradlance for non-zero albedo we use 

G(A*X*0) - — r (18) 

1 - r(X*Wj)A 

where A Is the ground reflectivity and 

G(O*X*0) - D(X,0) + S(X*0). . 

4. Analytic Representation of the Diffuse Actinic Spectral Flux 

Peterson (1976) has used the Dave Braslau multiple scattering code to 
calculate values, of the actinic flux as a function of altitude, wavelength* 
and solar zenith angle for an ozone thickness of 0.292 atm-cm, aerosol 
scattering depth of 0.204 at 302. 5nm and an aerosol absorption depth of 
0.034. He presents the values calculated for the ground level in two tables* 
one for zero albedo and the other for a best-estimate albedo. Starting with 
his zero albedo values* we subtracted off the corresponding direct actinic 
flux* D(X,0), to obtain the Inferred downward (since at ground level with 
zero albedo) diffuse actinic flux, S(X,0), at ground level. For the 
direct actinic flux Beer's law gives 

D(X,0) - H(X)e"'^®®‘^® 


( 19 ) 
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where t ■ ^ t.» where the t. are previously defined. Here we do not 
1 *■ 

differentiate the yj^'s since large solar zenith angles are not used. 


The downward diffuse actinic flux at ground level for zero albedo 
Is represented In a form completely analogous to that for the downward 
diffuse Irradlance (see Eq.(3)). The ratio 6) ■■ S (X,6)/i(X,0*) Is 
defined likewise In complete analogy with the Irradlance case (See Eq.(6)) 
as Is $ (see Eq. (7)). It will be convenient to designate the actinic 
case values of the GCS fit parameters and functions by piecing a tilde 
over the general symbol. 

Now, Instead of defining F In strict analogy with F In Eq. (8) we 
have simplified F a little to facilitate the fit to Peterson's data. The 
reason for this Is that the parameter g In Ea. (8) controls the variation 
of F with w^ and altitude and since we had only one value of w^ to work with 
and one altitude, the purpose of the parameter g was nullified. Thus, In 

** *** lu 

F we simply set g 0. M(X) Is defined In exact analogy with the Irradlance 
case (Eq« (13) )• 

The parameter values for this representation of the actinic flux are 
given in parenthesis In Table 1. Figure 3 displays graphs ofj^(X,6) 
andD(X,6) for selected values of X. One s^^.es a close similarity of the 
jS(X,0) with J(X ,0) in Fig. 1 in both shape and magnitude. TheD(X,0) 

•m 

andD(X,0) are quite different, however, because of the cos0 factor In Eq. (1). 

A comparison of the parameters In Table 1 Indicates that the y's do 
not change much: ckanges the most with a 35% decrease. F and F are 

both relatively small throughout the wavelength region of Interest. All 
In all, the normalized 0-dependence Is not much changed. 
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The greatest difference between Irradlance and actinic representations 
^ y (see Fig, 2), 

lies between /7(X) and A(A) j^We may first note the ratio " 7.95. 

The large relative value of causes hx) to fall more steeply as X decreases 
from 310nm than does/Y(X). Above 315nm the denominators o£ and/j[(X) 

are both nearly unity. Now turning our attention to the numerators* 

let us denote the Rayleigh and aerosol terms by R and A and take note of the 

following handy decomposition: 


<*•»** « 

R+A_R 1 .A 1 

- . + _ . ( 20 ) 

R + A R 1 + A/R A 1 + R/A 

■ mt m m. 

Now note that R/R ■ ^^l^^al " 2.19, A/A has the essentially constant value 
of .59, and R/A lies in the range 2.0— 3.5. The ratio' of the numerators is 
thus seen to lie in the range 1.6-1 ^8. Taking the effect of the denominator 
into account, one obtains that the ratio ^(X)/^(X) lies in the range 1.4-1. 7, 
and in fact the ratio /f(X)//((X) has a mean value of 1.62 over the wavelength 
interval 295nm-365nm (the wavelength region of interest) and median value 
of 1.64 there. 

Going back to the normalized 6-dependence, we observe another interesting 
contrast in the parameter values between the actinic and Irradlance cases, 
namely in the parameter A^^. The comparison of A^^ and A^^ is complicated 
by the fact that we set g "0. Noting that for X < 320nm, it is 

clear after a little algebra that Aq^ roughly corresponds to Aq 2(1 + Sq/w^) 
in this region. Noting that »T 2 for X > 340nm, it is then clear that 
Aq^ roughly corresponds to A^^ in this region. Thus one might expect the 
value of Aq^ chosen by the fit to be some sort of average of the values of 
the above two expressions which are, respectively, 23.34 and 3.285. This 
is Indeed the case, with «■ 
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5. Conclusions 

In the work sunmarlzed In F-ection 3 it has been shown to be possible 
to represent diffuse spectral irradlance In the 280-380nm region to a 
reasonable level of accuracy by an analytic model which displays the 
explicit dependences upon the ozone parameter w^ and the aerosol parameters 
T 2 and For a more precise representation of the Braslau Dave data, 
i»ii. to the 5% level of accuracy, the work of Green, Cross and Smith 
suggests that some of the parameters should be particularized to the 
specific Braslau Dave atmospheric model. It might be remarked that with 
the adjustment of only one aerosol parameter (^^ 2 * “ average GCS 

model (see Table I) successfully represented a large body of Irradlance 
data taken In the South Pacific by Baker and Smith (1980). 

The work presented In Section 4 is based upon a much more limited 
data set. It suggests that actinic flux can also be represented by the 
same equations. A detailed comparison of the ^(X,9) and.J^X,6) in 
Figures 1 and 3 Indicates that these two quantities are very close for all 
values of X and 6. Thus the difference between S(X,6) and S(X,6) is 
largely manifested in the differences between/t.(X) andy^(X). These 
quantities have essentially the same shape and the significant difference 
between them can be best expressed as a rough scaling factor of 1,62. 

Our work on the actinic flux is not yet completed. We have yet to 
test our model with other values of ozone thickness or at varying altitudes. 
Such work, which requires the generation of additional "data" using a 
radiative transfer code is now in progress. However, it is not un- 
reasonable to expect that the GCS dependencies upon w^ and y for spectral 
irradlance will also be efficacious for representing the dependence of 
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actinlc flux upon these same variables. As an Interim measure to char- 
acterize the actinic flux as a function of all Independent variables Includ- 
ing and y, it Is not unreasonable to usei(X» 0) -i(X.6). >l(X) - l,624f(X) 
and r(XyW^) ■ r(X,w^). The Interim and the final actinic formulae should 
be applicable in models of smog reaction, tropospheric chemistry and other 
Important photochemical problems. 
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Figure Captlonc 

Fig. 1. S(X,e) and P(X,6) for X - 297. 5» 307.5, and 330. Onm. The 
points represent the "data" for S(X,C) ; the solid curves 
are the CCS fit. The broken curves Illustrate 17(X,0) 
at the corresponding wavelengths. 

Fig. 2. M(X) and M(X). The points represent the "data" for M(X). 

Fig. 3. S(X,6) and P(X,6) for X ■ 297.5, 307.5, and 330. Onm. The 

•*0 

points represent the "data" for S(X,6); the solid curves are 
the fit arrived at In the present; work. The broken curves 
Illustrate P(X,0) at the corresponding wavelengths. 
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This book is a product of the collective work of a group of faculty 
■enbers of the University of Florida. Although e.ich writer Is a full- 
tlae teaching or researching neaber of a disciplinary departaent of the 
University, all participate regularly in interdisciplinary research con- 
ducted under the auspices of the Interdisciplinary Center for Aeronoeqr 
and (other) Atanspheric Sciences, known as ICAAS. Since 1970 ICAAS 
has conducted research projects pertaining to atnospherlc pollution and 
related phenonena. These studies are exeaqillfied by a nultldlsclpllnary 
aasessnent of an anticipated large Increase in the utilization of coal 
in Florida .which is now underway with Board of Regents support. 

Although nany ICAAS iwnbers ace physical scientists, its research 
has not been llEd.ted to hard science Issues. Rather. ICAAS enploys a 
broad gauge pidillc policy approach, bringing together the hard sciences. 
111 e and agricultural sciences, social sciences, econonics. medicine, 
law and other disciplines with the goal of seeking technical solutions 
to complex problems that are also socially and politically acceptable. 

ICAAS members foresee a period of severe social stress ahead as 
this nation is forced to switch to coal from oil and natural gas as its 
primary energy source. The forces are multiple, posing specific techni- 
cal. environmental, political, legal, and other kinds of issues that 
ere IndJ /idually Intractable and that in sum pose an extremely formidable 
barrier to the nation's well-being. ICAAS imBbers believe that a con- 
centrated nultidlsciplinary effort by academic scholars and researchers 
can help the nation get over that beerier. 

.Such efforts have many precedents. After the outbreak of World War 
11 ibe American academic coMunity rallied quickly to the nation's de- 
fense. In addition to specialized training programs, important rontri- 
butlons were made to technological developments, such as radar, missiles, 
cocket-aaslsted take-off aircraft, the proximity fuse, weather fore- 
casting techniques, operations analysis awthods. and nany others which 
helped win World War II. Beyond doubt, the atomic bari> is the best 
known World War II product of academic scientists. Although almost all 
of Japan's war potential ms destroyed before It ma used, the bomb 
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still is generally credited with avoiding the invasion of Japan which 
would have taken uany lives. 

Today, academicians, like all citi:;ens. are grateful that the 
nation is not engaged in a shooting war. Nevertheless, the Independence 
of the nation is threatened so severely by the Impending energy short- 
fall that President Carter has called it the moral equivalent of war. 

The relatively slow development of this "moral war." which truly is an 
economic struggle for survival, affords the academic co—«lty with yet 
another opportunity to serve the nation. The nature of this war and 
the general direction that the United States must now take have be^ 
outlined in Project Independence (1974). The National Energy Plan (1977). 
and the Camp David plans (1979). Although wide areas of disagreement 
about what should be done still remain, both RepubUcan and Democratic 
administrations have agreed that the essential alms set forth in these 
plans must be achieved: (1) import less foreign oil; (2) conserve ener- 

gy; (3) to live within the nation’s resources, make coal the major al- 
ternative to foreign oil. 

How can the academic co-unity help achieve these goals? In gener- 
al. it is not as well-equipped as large industrial and governmental lab- 
oratories to handle major engineering developments, but it can originate 
new ideas, evaluate and extend innovative technology, and carry out early 
phases of develrpment work. Perhaps most Important, a University facul 
ty. being a repository of experts and scholars of many disciplines, is 
unusually well-qualified to perform multidisciplinary assessments. 

This book singles out one topic, namely Coal Burning Issues, as the 
focus of its attention. The writers have examined the technical, medi- 
cal. environmental, legal, economic and public policy issue, that must 
be addressed as the nation increases its dependence on coal. The book 
does not purport to supply answers; too much .»rk remains to be done 
such as a more in-depth examination of the apparently most feasible solu- 
tions. Instead, it is written with the hope of accelerating examination 
of a series of critical, long-term strategic and short-term tactical op- 
tions. This must be done without undue delay because passing time Irre- 
vocably closes options- 
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This book employa no single style. In pert it employs qo— clcacivm 
«»lyses «>d in psrt it examine. •oci»U poUtiesl sod legs! choice. Im- 
fluencing the nation*, energy future. Bo attsivl 

to eliminate all dl— greementa. For example, the reader rntj mote di- 
vergencies among so— of the chapters in estl— tes of cort and other re- 
sources. These esti— tes have traditio— Uy varied by large factors de- 
pending upon the perspective, of the e.tl— tor. and diverge— la. are 
char-terlstic of the technical pshUcatlc- examined In the cour- of 
this study. Accordingly, they do —t detr-t fr— but Indeed enh— ce 
this work, which al- to hi^ight both -ttled and disputed f-t. and 
viewpoints that wlU Influe— e Inmlncnt public poUcy choices. 

In the past. A-rlca*. «rademlc coumunlty has played a vital role 
in protecting the — tlon in shooting — rs. It see— certain that — 
less a role is to be played in wl-lmg the -onomlc -r of wirvival that 
now entangles the -tion*. future. The writer, of this bock hope that 
the pre-ntation — de here point, the —v to — lutlo— to — of the 
coal burning l.s-s. The next goal of the group i. to help supply — 
of the detailed solutio— to the energy, environment end — o— -qr 
proble— related to the 1— reased utilisation of coal. 
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CHAPTER 1 


INTRODUCTION AND SUHHARY 

interdisciplinary Center for Aeroncsv and (other) Atmospheric Sciences 


1. IHTRDDUCTIOM 

n . rf tt . oU "•'<“1 «" »■ 

.utiT a» 1- ^ 

tb. ~U. ™. 40 »7 >■“ «• “ 

■sde and old -eUs are rejuvenated -ith ne- recovery -ethods, no kno-1- 
edgeable per-n. d««y that ulti-ate exhaustion of stored fossil fuel, 
is Inevitable if consuuption continues at present rates. The urgency 
the proble- -as underscored In 1979 hy three events that have had a 
«Jor inpact in bringing the energy supply situation to the crisis 
stage. First, the revolution In Iran -^.planted a govemnent friendly 
to the united State,: and that ha. recently -ide extensive purchases In 
return for oil -ItL one that is -ragonistlc to this nation. Secon . 
the Three-Mile Inland nuclear reac .. accident has had a 
on the grouth in n«her of nuclear rc tor electric 7 '' 
thus undercutting a -aln alternative to the ue* of oil «,d gas ttat the 
nation had been counting on. PubUc concerns over storage of radlo- 
«:tlve uastes the proliferation of nuclear technology have no- been 
aug-ented by concern a. to the safety of current reactor, and t^ 
adequacy of training provided reactor operators. Third. OPEC o p 
increase, have forced the United State, into virtual econo-lc -arfare 
-ith -any of the oil producing co«itriea; a -ar that the nation 
appears to be losing, if the rising price of gold, the rising interest 
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rate, and other signs of Inflation any be taken as indicators. 

Ihe urgency of the need to reduce purchases of forei^i oil and, 
therdjy, help stea the net outflow of dolUrs to pay for iaported 
energy is now so acute that alaost everybody acknowledges it. 

What can the citizens of the United States do? The aost obvious 
ttdng. albeit a little old-fashioned, is to take stock of our natural 
resources. particuUrly our renewable resources, and use Aaerlren com- 
mon sense to find a way to live within the nation's aeans and re- 
establish national self-respect. Highest prlorxty nist be given to 
conservation and the use of renewable resources since these actions do 
not Involve any drain upon "energy savings." Nevertheless, phasing in 
a renewable resource node of life wiU take tine and could be rather 
painful. Fortunately, the United States has a goodly supply of coal 
that is estlMted to awunt to 302 of the world supply, enough to sus- 
tain the nation for a hundred years or so before a new energy resource 
is discovered. Accordingly, coal is a resource whose use aust now be 
pronoted. On the other hand, coal poses dangers of polluted air, dirty 
water, aangled earth, crushed bodies and bUckened limgs. t— g .. „f 
these hazards were created a generation or two ago before oU and gas 
becane so abundant and cheap that they ousted coal as the fuel of 
choice in Antrica. The nation is now co^ieUed to return to coal and 
find out how to utilize it without resurrecting those bad laages. 

Modem technology is the key. 

This book exaiaines the potential role of coal in the United States 
in the coKlng decades. It looks at coal utUlzatlon broadly to 
burning and other direct uses. Uquefaction and gasification. To assure 
conprehension. the perspectives of scientists, engineers. systeM 
analysts, nedlcal scientists. Uwyers. econoiasts and other experts 
have been eiqiloyed. ParticuUr attention is given the interln role 
that coal Bust assume to neet Aiserica’s energy needs in the next 
generation or two. This will give a breathing space for science to 
ascertain whether or not unllnlted energy resources such as fusion 
■achines or breeder fission reactors are feasible, and for society to 
decide whether or not they are acceptable. If afflmatlve answers to 
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both questions are not forthcoalng after a few decades, society nast 
then use its renalnlng coal and other fosail fuels to to live 

entirely on soUr energy la direct or iadiract fotM. 

TWa book has been writtaa with the asawaptloa that society does 
not know what the future holds with respect to these loag-rsnat alter- 
natives. This will be disconfortlag to sow readers who have flrwr 
visions about the future than do tba scientists and other experts who 
wrote this book. This underlined the Inportanee of the coUective 
point of view shout the use of coal expressed by these writers-^osl 
buys tlue for this country and the world until the ^pe of the long- 
range future becones clearer. Thus, this book coneeins itself with 
the best way to utilize coal fro. the viewpoint of our energy needs, 
our envlronaeatal safety and the health of our econo^. 

II. QUAHniATIVE FACTS ABOUT BIEIGT SUPPtlES A«> COMSOIVTION BATES 
The old saw about history's contlnuaUy repeating n— if has 
obvious validity so far as the use of coal is concerned. The curves 
of Figure 1 depict how during this past century or so the U.S. has 
shifted frou a renewable resource, that is. fuel wood, to coal, and then 
to petrolew and natural gas. Each shift has been spurred by factors 
such as availability, convenience, and econonlc advantage. Since coal 
is stlU abundantly available it obviously ..st have sow reUtlve 
dis.dv«.t.ges that wst be coped with. Alwst everyone has heard rtnuit 
Pittsburgh's reputation as "the awky city." a nicknaw that no longer 
appUea. and the descriptions about the coal swke pollution of 
are even «,re fear«*e. Perhaps fear, both legltlwte nmd apocryphal, 
about the bad aide effects of burning coal even deli^ed the onset of 
the Industrial revolution in England. History recounts that in the 
year 1306 King Edwari 1 of England feared it so aucfa thst he decreed 
punlsfawoc of death to people who burned coal. Ultiwtely. however, 
the power of coal was not to be denied ev«> in England because history 
has recorded that about four centuries Uter that nation used coal to 
fuel the industrial revolution that wde it a great world power. 
Sinllarry. coal has been a wjor source of energy used la the United 
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are charted ia Table 1. . -i _ 

The world Is now on the threshold of shifting hack to coal as 

«ijor source of energr- Convenience and environ.ent.1 adv«.wge no 

longer are the driving forces; this ti-e it in one of necessity ^ 

ueirknown to everyone, petroleu. and natural gas are 

. •• a ualor source of energy sowetxwc 

will by all accounts he depleted as a waj 

rj.i». -u .. .h. -h 

drl... «. U. .< — ... ».l»i h»... “* 
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driven downward. _ ^ 

Although this book is about coal use. the writers do 

leave the erroneous Iwpression that energy conservation is uni-portan . 

It trfces no *,re th«. a sU|^it digression to de*«strate that conser- 
vation can offer substantial relief that, when coupled with Increase 
coal use. can help return the United States to n status of energy 
independence and sustain it there. Per capita ener^ „ 

the heated States is so high by any standard that the very fact s^st. 
that censervation -ay he of «Jor guantitative Inportance. For exa-ple. 
ss shown la Figure 3 (McPherson. 1965). Anericans far exceeded a^ orter 
of the world-s people in hoih capita energy consu-ption ^ 1«^ 

in 1965. Co-parative data for 1974. pictured in Figure 4. show tha 
in IVOJ en-rev in the United States again 

although the per capita consunption of energy 

exceeded that of any other nation, the equivalent per capita inco-e in 
«ay IndustrlaUxed nations had overtaken or even exceeded that en- 
joyed here. Taken together, these factors suggest that * 

to A-erlca wight produce substantial energy savings without curtail 

the standard of living. 
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Table 1: Som» Biatorioat Coal Burning Landmjrks 


Tine 

Country 

Person/People 

Use or Context 

ISOO BC 

Halea 


coal used in funeral pyres 

1100 BC 

China 


heat for howes 

950 BC 

Israel 

King Solowon 

mentioned in Bible 

300 BC 

Greece 

Aristotle 

noted disagreeable smell 

121 AO 

England 

Hadrian 

used In dwellings near wall 

852 AD 

England 

Anglo-Saxon 

coal used in rmtal payment 

1306 

England 

King Edward 

decrees use of coal punish- 
able by death 

1200 

North Asttrica 

Pueblo 

coal used in pottery pro- 
duction 

1679 

North Aaerica 

Fra Hennepin 

observed black mineral 
along Illinois river 

1694 

England 

Clayton 

first production of coal gas 

1750 

England 


coal fuels industrial revo- 
lution and rise of England 

1770 

United Colonies 

Washington 

noted coal mine on Ohio 
river 

1771 

Pennsylvania 


anthracite discovered In 
eastern Pennsylvania 

1840 

U.S.A. 


coal Industry grows 
1,000,000 tons mined 

1850 

U.S.A. 

Gesner 

real liquefaction 

1890 

U.S.A. 


electric steam generators 
expand use of coal 

1940 

Germany 


large-scale liquefaction of 
coal by Berglns process 

1950 

U.S.A. 


Middle East, Venezuela oil 
flood world markets 

1973 

U.S.A. 


oil crisis begins 

1974 

U.S.A. 

Nixon* 

Project Independence 

1977 

U.S.A. 

Carter* 

National Energy Plan 

1979 

U.S.A. 

Carter* 

oil crises worsen - Camp 
David proposals 


*propojes aajor switch to coal utilizatloa 
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Ihl* My be illustrated by co^wrlng the fuel efficiency of 
AMricaa cars to that of inported cars. In 1975 the average doMstic 
car used shout five gallous of gas to travel the saw Jis tafw as the 
average inported car uould go on three gallons. Hence, auch petroleia 
could he saved hy nerely converting to fuel-saving cars. This will 
occur if present federal prograna successfully bring the average rate of 
fuel use of the doMstic fleet down to that of inported cars. The waste 
transportation apparently has its industrial counterpart if 
conpariaona between Aaerican and Nest Gernnn industrial fuel consunption 
data are a reliable guide. For r::.aple. the West Gernans use about 
20H less energy per unit of prinary netala produced than does Anerican 
industry (SM, 1975). Other najor energy intensive industries of West 
Cetnnny have even greater relative efficiencies. Clearly, a role for 
better industrial technology in energy conservation is evident. Hhat 
happens in industry and transportation seens ti represent <iolte faitb- 
fuUy the relative energy voracity of the AMrican way of life. The 
***** *•■■■* energy Inpoaed by the average Anerican coapared to that 
for people froa other countries has been noted by Sail (1979) who coa- 
pares the per capita aaonnts of energy consuaed aeasured in Ulocalories 
(heal per day) of the developing world with that of the USA. 
tnia work shows that Aaerican per capita energy consuaption is exorbl- 
!■**'*• About 215.000 kcal per day whereas nany of the world's 
developing countries get along on about 10.000 kcal per day per capita 
with little fossil fuel a u gn en tation to solar energy input. In China, 
for esaapls. the energy use per capita is approninately 13.500 kcal 
^***A» Anerlcans could help alleviate the energy crisis by 
■owing the per capita energy use dowward. possibly to 100.000 kcal per 
day. the approninete energy use of other western Industrialized nations. 
If. for ezaaple. we could Just lower the rate cf growth of energy use 
froa 3.5Z to 2.3* per annna we would save 20 ailllooa barrels of oil 
per day by the year 2000. 

afcjactlve of this book, however, is not to dwell upon the 
energy conservation aeasures that the writers luianlnnusly subscribe to, 
but rather is to ezaalae phenoaena associated with a transition to coal 
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use. particularly in the United States. The use of coal in this 
country any be what will aake it possible for Aaerlcans to aalntaln a 
high standard of living, leaving the foreign oU to alleviate the 
plight of poorer people In the world. 

Fortunately, the world and especially the United States have large 
supplies of coal to consiae. In gross dlaensions. the world's cial 
reserves are Ulustrated by the data In Table 2, extracted froa Peters 
and Schilling's (1978) appraisal of worU coal resources. The data are 
presented in aetric tons with the energy equivalence of hard coal (HC) 
rather than in gross weight. Hard coal here encoapasses anthracite and 
bitualnooa and brown coal (BC) includes sub-bituainous and lignite. 
Gross supply data also have been adjusted to reflect only supplies that 
are econonically and teclmlcally recoverable. Hence, other estinates 
of coal resources could vary substantially froa those shown. Be that as 
it any. this table indicates that the United States has aore recoverable 
coal resources chan any other nation. The last two rows give the U.S. 
production and export cates in 1975 in wUllon tons. 

Table: 2. EaonanoaZZy and TeahnioaZZy ReaooerabZe 
Coat y-sotepoae of Leading CoaZ Countriea 
in Gigatona {Vfl metvio tone) 




HC 

BC 

Total 

Prod 

Exp 


Horld 

492 

144 

636 

2593 

199 

1 

USA 

113 

64 

177 

581 

60 

2 

USSR 

83 

27 

110 

614 

26 

3 

China 

99 

— 

— 

349 

3 

4 

Great Britain 

45 

— • 

45 

129 

2 

5 

West Geraany 

24 

10 

34 

126 

23 

4 

India 

33 

— 

33 

73 

■ 

7 

Australia 

18 

9 

27 

69 

29 

8 

South Africa 

27 

0 

27 

69 

3 

9 

Poland 

20 

1 

21 

181 

39 

10 

Canada 

9 

1 

10 

23 

12 
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Peters and Schilling also have predicted how rapidly coal will be 
put Into production in the coning years. As shown in Figure 5. the 
United States is expected to lead the world in this process. Whether 
or not It does will depend upon whether or not and how the Issues 
exanined In this book are resolved. In any event, it seens certain that 
a substantial shift-over to coal is coning. Already this is nenlfested 
by the fact that virtually all new electric power plants purchased In 
the U.S. are being designed for coal rather than oil or gas. 

The fact that shifting to coal opens up the world's last known 
najor resource of non-renewable energy to consunption. with an appetite 
that will Inexorably lead to ultinate exhaustion poses a noral and re- 
source issue for bunankind that has yet to be squarely faced. What will 
be the plight of those unborn generations that cone along when the re- 
source is depleted? What should we do about the underdeveloped nations 
who are H—inding the right to develop? This necessarily inplies novlng 
away fron their current equUlbriun with renewable energy (Snll. 1979) 



fig. 5. Survey of the Future Trend in Production 


toward a dependence upon stored energy. Certaliay if China's bUlion 
people were to nove up to Anerlca's per capita energy eonsnaption rate 
the world's stored energy resources would be wiped out In short order. 

This world-wide social dilen* has yet to surface in fnU force. Unless 
a new, ineAaustible energy resource is discovered, considemtion of 
equity suggests that sons sort of baUnce »st be struck in which Chose 
with nor*, use less and those with less, use nore. Boiiever. the aUo- 
cations to each nation of energy resources, property, population Units 
and enlssiott Units to control adverse cllnete changes miaes conplex 
international issues considerably beyond the sc. pe of this work. 

Another issue whose conplexity nakes it beyond the scope of the 
present study is the question of nonopolistlc practices la the energy 
field. Such practices by the OTK «U cartel obviously have contributed 
to our energy crisis. Voices are now being raised that najor Anerlcan 
.-n-peniM hsve established an energy nonopoly pamUellng the OPEC car- 
tel. However, it is an open question as to whether this cry identifies 
a ■■ ■■I problen or alnply diverts the nation's attention fron fwndsMntal 
problens: the inevitable eidiaustion of the O.S. oil and gas reserves 

and the need to use coal to save our econo^r. 

The existing coal suppUea of the USA are large enough, even if 
used up at a voracious rate, to provide tine for the developnant of even 
nore plenteous energy resources that nay exist, such as nuclear fusion 
or breeder fission, or to learn hew to live on renewable energy re- 
sources better than the developing countries. Beferring again to 
Figure 5. the reader can sec that the actual 1979 coal protection is 700 
■illion tons per year in the United States. At that rate, the eatinnted 
Anerlcan reserve of 177 Gigatons of econonicaUy and cechnlcelly recover- 
able coal would last 250 years, and even if the eonsMptian rate were to 
triple, the supply would last nearly 100 years. This Mderacores again 
the central purpose of this book; nanely. to exanlne the isaues and 
problens that nust be resolved to sake coal a sags and eeviroansntaUy 
acceptable energy workhorse for the USA during the coning 50 year or so. 
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III. OIKIGI OMITS AND THE DIMBISIOMS OF THE FEOBLEM 

TIm units us«d to describe energr qusntltles eery grestly. even 
emg professions! people. ftMsn energy connonly is neasured in con- 
s^tion of kllocslories (keel) per capita per day. For exanple. under 
■oderate exertion the average Anerican sale uses 3.220 kcal per day and 
the average Anerican fennle 2.320 kcal per day. By contrast, heat engi- 
neers describe energy In terns of British Thecnal Units (BTD). physi- 
cists use the kilojoules (kilouatt-sec) . electrical engineers use kllo- 
untt-hours. petroleun engineers use barrels of oil, gas engineers use 

of cidtic feet, md nlnlng engineers use tons of coal, either in 
ahort-tons (2,000 pomds) or netric tons (1.000 kilograns). All de- 
scribe the saM basic qnsntlty. but In terns that have no readily ap- 
parent conversion equivalences, such as the trivial conversion of, say, 
feet to Inches. This conwnlcatlon hitch has been nade even worse by 
the growth in use of a nunber of new words to describe energy resources 
in global terns. Thus, the quad - 10^ BTU, the Gigaton - 10* tons of 
coal equivalent (GIC). have cone Into use and now conpete with older 
terns sneh as -ilHnna of barrels of oU (ISO) and trillions of cubic 
feet of natural gas (TCF) to describe energy resources. 

To supply a connan base for cnnmsil rating about energy quantities, 
the InfotMtion la Tbbles 3 through 5 is provided. Thble 3 nerely de- 
fines the wilts; Table A displays aone of then la a fashion that nakes 
conversion anong nicro-units relatively easy; and Table 5 displays 
others in a fashion that does the sane for nacro-units. The BTO is used 
to tie Tables 4 and S together. 

Am a gaMirlretlve o v e rv i ew of the dinenslons of the United States 
problsn In terns of these units. Table 6 gives a detaUed picture of the 
H.S. energy supply and desnnd fron 1960 to 1978. Headers nay find it 
helpful to nark these tables for future reference. Right n<nr this 
nation's coni prodnctlon capsbllity far exceeds the current denand and 
production could go to as U.^ as 2 Cigatons by the year 2000. Thus, 

could replace oU and natural gas as our naln energy source. 
Uhether or not it should or would is the nela concern of this book. 
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The foUowing sectloa of this chapter conalats of salaries of ' 

Chapters 2 to 18 of this book. Chapters 2 to 8 lartcly consist of \ 

descriptions of our resources, present and future coal tnchnologies and 
current coal issues. Chapters 9 to 14 describe knonn or possible 
environnental probleas. FinsUy. Chsptets 15 to 18 disenss and 
describe possible asses a nen t techniques ubich ni^t be used in public 
policy decisions related to coal burning iss u es. 

In addition to the references cited in this chapter ue list below 
a brief bibliography representing uety re!*.enc norks on energy in genetnl 
and on coal utilization in partlenlar. 


baeatlM Offtca at tlw rRsMaac, Wt, Sttn 
Cun— r Prlaclae Offle*. UMUacua. •• C. 


r Mlcjr mi n — l ae. *• 


Sen— m. W. IMS, Egamlc_ ^li|ai»t ff 

Chapter t. pp. lam State Balraralty Piem. mm. Ima 

g”. **j **pmlml ef UmU Omi I 


■eiteiiite. B.. U». m letteiecti m te Cml IhchmlmT. maiaaic 

Cram. A. T. at el.. »», the Street See ef Caal. Praapecta aai 1 
r Iji - ef the eaitai State*. Office ef Tedmlact iaaeaam 

Beela. R. B. et el.. U79. «vi ■tl«« y3gjB B egBSy i 

...tea tm the Soath. Oak Uipe Batiaaal Lat-acaCarlm Bapart 

Bepartaaac ef lecxer. B. S. M7P. Ihel Bm ict Tl 1 til F 

Battel. B. C. aai Boaati. J. B.. Iftl. Bm Beaiap T-e— 


I B. B. et el.. 1919. reateet the Seat theacr 1— 

■dllaear PhhUahlae Caapaap, CaabcliBe. Bam., pp. 21. 
Uaaeci. J. B., 1S7P. Coal, the Berli Brah thac Baek. meU 
lU.. pp. SM-Stt. 


teftaam, B. 1.. U7B,- 


k. tm SMCiaai BelaheU C. Bm ' 


HcBaal. B. B. aai Blelaaa. C. P. et at.. U11. Baaatam Caal 1 
cadme, Bm Tack 

Btakaeih. B. aai tereia, et el.. 197P, Baetkl mtme. Bajertj 
sawel. Baaim Beam PiAliahlnt. Bm tack. 
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Chaptar 2. Coal Aira H abt ll tv and Coal Htnlaa. by Ohaalan and 
Fardablahah. ataanirlzaa tne U.S. coal resource base by geographic lo- 
cation, access (l.e. . deep or shallow nines) , and clasatficacion. A 
ganeral asaesanent of denand projections is based on the assunptlan 
that elaetrlcal pow e r genacatlon be the priaaty user of coal over the 
next two to three decades. 

Coal is the nost abandant non-renewable energy re&'onrce in the 
United States which has approalnately 28Z oi the world's total coal 
anppty. U.S. coal, which represents 70Z o£ the non-renewri>le energy re- 
aoarce la the coantry is located In three distinct regions: Eastern, 

Central and Weatem. Bitanlaons coal conprlses A3Z of the resource 
base, sab-bitanlnons 27Z. lignite 28Z and anthracite IZ. The sulfur 
content of Vestem coal is lowest (generally IZ or less) with Central 
coal higbaat at 2-AZ and Eastern coal at nost 2Z or less. However, the 
heat value of coal per unit weight is. in general, lower for the 
Westein coal reserves. 

A mmmmry diacueslon of nodara nlnlng tacholques and the factors 
that give rise to the use of a particular nining nethodology are also 
preaentad. The role played by tacfauologlcal developnent throu^ the 
evolution of nee bine ry and agulpMnt used in the nlnlng syaten is dis- 
casaad. lha anvironnantal isanes arising fron coal nining are exanlned. 
Since coal crashing, washing, drying, ate. is often required, it is 
con al d ar ed an integral part of the nlnlng, and coal preparation is also 
covnmd briefly. Siaea safety plays a algal f leant role, the problens 
mlatad to safety ace dincaa a ad. A related topic which is of paranount 
inportanca. sanely the o c c u p a tional hasacds associated with coal nining. 
is ravl awad. It la saan that since the anactnant of the 1969 Federal 
Coal Htna H e al th and Safety Act. nine fatalities have been redact'd 
appreciably. Howavar. no gtantly naaaarabla rednetlon In the rate of 
dlnahling Injariea has ocenmd. Pigaras ace tabulated which show 
nlners ate likely to suffer fron occupational health and safety hasacds 
at a cate higher than other occnpatlons. 
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Chapter 3, An Enernetica Analysis of Coal Quality , by Alexander. 
Swaney, Rognstad and Hutchinson uses energy circuit language, a netliod of 
•y*t*wo analysis, to detemlne the quality of energy enbodled in coal. 

Three distinctly different approaches ace ei^loyed to calculate coal's 
energy quall^. 

The first aethod enployed to detemlne the energy eabodled in coal 
uses a theoretical nodel of the coal fomatlon process- As a startlns 
point a detailed description of the fomatlon of tHe various ranks of coal 
is given. Fron photosynthesis, atnospheric CO^ is transfomed to plant 

to peat via niccobial deconposition In water. Then through 
catagenesis (autoclaving geological processes), peat is transfomed to lig- 
nites to bltunlnous and anthricite focus. The geological and ecological en- 
eegies Involved in these tranfomatlons are evaluated. 

The second nethod used to detemlne coal quality conpares the energet- 
ics of coal burned as an electrical generating fuel to the use of wood and 
oil In electrical pewer production. These two fuels were chosen as a ba- 
sis for conparlson because the energy value of wood is easily traced to 
the energy flows of natural systens and the energy in oil Is the largest 
oonponent of the contenporacy fossil fuel energy base. 

The third approach to estlnatlng the energy quality of coal is based 
on the energetics of converting coal to higher quality synthetic gas. It 
is postulated that this will give an upper range to the energy quality of 
coal, while conversion to electricity wUl appcoxlnate the lower Unit. 

It appears that 1.36 — 2.2 themal calories of coal are required to 
do the sane anouat of work that one calorie of oU or gas can 
In effect, this aeans that the nagnltude of the world's coal reserves ««»«~. 
expressed in themal units, should be reduced by 26 - 55Z. 
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Chapter k. Coal Tranaportatlott . by •• Capehart, O'Coimor aod Deaty 
Che li^ortant H"fa« In the overall processes of supplying coal 
to latiatte users. Coal Is shipped by a variety of nodes and nay even re- 
quire nultlple nodes, were 65X by raU. UX by water. 12X by truck and IX 
by other nodes (slurry pipelines, conveyors, etc.) Eleven percent was 
used at nouth plants. Future expansion of the coal transportation 
network is clouded I •- v ssent federal policy which appears to eaphasise 
short-tem expansion o.. ns'-ural gas use. rather than coal. If the coal 
transportation systen is ;o be expanded, the rail systen which can nost 
easily be i»» panH«wi nust carry the najor burden of growth since it is the 
nost easily expanded transport systen. Expansion of the coal transporta- 
tion network will result in najor environnental and social problens atten- 
dant on both construction and operation. Air. water and noise pollution 
acconpany the different nodes in varying degrrs. Social inpacts such as 
coanunlty disruption and safety also occur. 

Operation and future construction of transportation nodes ate greatly 
affected by federal regulations and reguUtory boards. The extent of reg- 
..la Hnna concerning the different nodes varies fr«e alnost no control to. 

In sone cases, a conplete roadblock to expansion. Brpansion of coal 
transportation will require the ex a nin ation of najor public policy deci- 
sions. Environnental. social and econonic problens ate all Interrelated, 
and nost be considered In a systens approach. Conplex cost/benefit studies 
will be necessary to help deteialne these public policies. Finally, eco- 
nonic conparlsons of the various nodes of transportation are needed in or- 
der to nake optlnun econonic decisions and expansion plans. This is also 
a conplex area since econonic decisions are not independent of environ- 
■entalg soclftl uid regnl^tory rdiolTCMBCSa 

The future of the coal transportation systen la still very uncer- 
tain. Until questions of future denand are answered note conpletely. 
few mjor Investnents in expansion are likely tc be wade. 
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Chapter 5, Coal Enmina Tachnoloar. by Schnaidar and Emm la de- 
voted to coabustian ae p aets of the diract nan of coal. Ihn diacnoalan 
centers prinarlly on coneentratad nan aneh aa hy olacteic n tillt lao. or 
Urge indoatries and dllnta use. as U roa l den ti al coal hnralng uhleh 
probably wiU bo an envlroonantal preblan U popn l a t od ragftiiaa. Conean- 
trated tochniqc.ua provide for mesmamf of aeaU for the ronoaal of baxnfnl 
pollutanta. The Inportant paranatars la coal coabnatUn am the baat of 
eoabnation. tba burning t enpe r a t a ra . the o xid a n t control, homing tiaaa 
and tba coal qnaUty. Ihesa paxanatom control the e nerg y ontpnt aa uaU 
as waste renoval state and nasoons by-product production. 

Historically, tba first typo of coal bnmar was the "flood bod" 
burner wbera fixed bod refam to tba location of the coal burning sane. 

Here heat riaes to an ov ecb aad boiUr to produce staon. In tbaaa In- 
stallations the stoker for feeding and coal aad ash renewal involveo the 
najority technology. 

Snspenalon bumam in idiicb coal dust is injeetad into a fim bm 
surrounded by boUer tubes r apraa e n tad tba next Inprevanant in coei 
burning techniques. Inprovenents of Injaction tacbaology led to tba 
present day "cyclone" furnace. Snepeasion burners have a very high tenp- 
eratnre which offers good baat transfer aad high off iciancy. Oafoctmate- 
ly. tbqr produce large anonats of ge a a n n a HO^. 

"Flnidlxed bed" bumem derive their nans fron the llqeid-like prop- 
erties of tba coekestlng naterinl. Hitb Inproved heat transfer techniques, 
fluidised bed burners can be opemtod at lower tanpe r e t e r ae thno pro due ■ 
Ing nuch less MO^. yet producing tba aana anenat of staan as a bi^ 
tenpemtum suspension bnmar. fdrtbemom. chantral tachniqnao in tba 
fluidised bed can be used to renowa SO^. Otba; ta ch n iq pa o for the di- 
rect use of coal which am still in the developaant stage in cln da n aga et o- 
hydrodyaanic (MID) ra a c tom and gna tarbine-staan tnebine hybrid ayatana. 
In the latter the gas tnrblna conld be snpplied by an MHD nossln or water- 
gas coabnstor. The gas turbine nduMst ■rvnld than feed a convontional 
boUer to produce etean. Ibe established techniques and eipe rl nen t a l 
teclaiiqnes are aU described aad jse e eed in tbs chapter on coal becn- 
Ing technology. 
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by Htrer .nH I^e. 1 . devoted 

to eoel ccverslo., the pcoceMo. of «ki„g clean burning liquid and 

fual. fro. coal. Tb. baalc concept of coal co.ver.lon chenl.try 
i^anrlch t« coal ulth hydrogen In ao« for. to produce l^drocarbon 

1th rlatlvly large hydro*, content. Thl. can .e acco*,ll.h.d 
by th. raactlo. of col 1th hydro*, or ate«i. or by re*„rlng e«e.. 
carbon fro. col to give a hydrogen-rich product and a carbonaceou. by- 
product. The ll*id fuel, produced in the llquefactl.^ proce.. range 
Cro. light gaaoUne to hMvy tar.. Uquefactlon and gaaiflcatlon reac- 
tor. involve ««y ..^rlapplng precede, and .any reactor, produce both 
ga*. and ll'juld product.. 

Seven co^erclally operating gaalflcatlon proce..e. are nou in 

ctcge- Thews gaalflcatlon pro- 

««.. can ba categorin In ter* of rwustor type and product. for*d. 

^ ra. *nthe.la ga. produced In «.t gaalfler. can be u*d either a. a 
fn or a. a feed atock to nek. *nthetlc ..tun gaa. wsthwsol. hydro- 

•"* ** hydrocarbon fuel., via the Fisher Tropwsh ayntb- 

ala. 

Col cow«r.lo. to pathetic fual ha. the grant potential to acco.- 
*>d.t. cleaning up proca*.. and could prolde an l*«rt.nt part of the 
clean gaMou. and Uqld enargl.. of fuel, needed for the next generation 
or tuo. Iha tachnlogy 1 . w.Unbl. for prodnclng both g.^ ..d 11- 
goM fael. fro. coal and only aconole conldaratlona. tl. capin in- 
ta*iv. nature of col cowr.no. plant. ..d the uncertainty a. to the 
CMC earn controlled price of fon*. oU hu Inhibited ow^ l«rge-*1e 
^•••lop*nt of coal eonvaralon In *i*rlra 

Tha poolhuttlan of converting oil ahalas. tar -a- m heavy pe- 
roeldnal. Into Uqld or ga*ou. .^thatlc fuel. .1* cowitltute. 
an Inportant lean., ih. raUtlvo acowsic. of the* *1.1. va. col- 
fcrlvl pylnal. iu datarlna th. pyathatlc ful. uhleh have the great- 
-t ah ort-rang. nmi long-rwig. prol* m the United State.. The envl- 
w—can Inpnct of tha axp«.d1 col ming and dlopon of large anoont. 

*** prodod In col eonvaralon la another la*>rtant n r 
related to coi convanion. 
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Chapter 7. Tyhnological Innovation., by Hanson et al . simrlzes 
«).. technological advance, which Ight anellorate. postpone, or even 
ellnlnate the energy shortage. 

Integrated UUUty Sy.te* (lUS) a. a concept originally Involved 
the tactical use. as dictated by ecoodc clrcuwitance. of gaa. oU. 
coal md «,Ud waste. However, because of the cost and unavallrisiUty 
of petro-fuela, the IU3 concept has evolved Into one In which coal U 
augaented by processed soUd waste. Never tiseless. such units have cer- 
tain advantages. Including financial, which can wake the. the dt of 
choice. 

E lectrically Propelled Autounhfle. « old technology, perlts the 
energy burden of travel to be shifted fro. gasoline to coal. Batteries, 
which have been greatly l^sroved. are charged by coal-produced electri- 
city thus freeing the car owner fro. the dependence on petro-chdcals. 

Further, by charging at Ight. the load on the power plant becowis 
better distributed. 

Ur Pollution Control Technology cold sake extensive burning of 
coal acceptable to society by proper control of effluents and residues. 
There Is continuing research In all phases of the tectoology of pollu- 
tion control. PartlcuUte. and SO^ have been given greatest attention, 
but HO^ control is also under developnent. 

Col Clealng nakes It possible to ellnlnate noxious or undesirable 
by-product, before the coal 1 . actually burned. Techniques have been 
developed for rewrrlng Inco^ustlble ash. pyrltlc sifur. and water fro. 
unprocessed coal. 

Off -Shore Fewer Plants orlgtnaUy proposed for nuclear power plants 
Ught be applied to col-flred units to Isolate and buffer the power 
faculty. Since al«.t half of today's power denand exists within a 
200 nUe strip long our various coasts, the concept Ight find ex- 
t^sive applications. 

Co aj ^ Plant Siting Technlyn.. in general provide the opportunities 
of cm plants to naxlnlxe operational efficiency, nlnllxe electric 
trananlsslon losses and environnental Inpacts. 
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Chapter 8. toter Reaonrces, by HelvUle and Bolch, sogsesta that 
two of the prlMry cooatralnta on the utilization of coal energy will 
be water availability and the potential for pollution of water tlMt 
would othenrtae be available for other uaea. «n analyal. of the avaU- 
able freshwater reveala a breakdown of 95X grounAmter. 3.5X lakea. 
swaupa, reaervolra. and river ehannela. and 1.5X aoUd nolature. The 
groundwater atreaaea are enphaalzed and are ahown to be the noat 
threatening. 

Surface water poUution la eaaUy detectrfile. The Mat nenaelng 
characteriatic of grouidwater pollution la that damge la probably Ir- 
reversible or at least corrective ttae scales could be ceasured In tema 
of years or aven generations. 

Suggested actions for reduction of water resource constraints are 
listed. 

Ale of the key Issues In the Increased use of coal will be the al- 
location of water resources to conpetltlve users. The conpeting water 
denands are enuaerated and discussed. The differences between point- 
source and non-point source poUutlon are discussed as they apply to 
coal burning facilities. 

The possibility of utilizing to good adv«,t.ge the water content of 
lignite, which is usually a problen, to faciUtate coal conversion is 
briefly discussed. 
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, Aa s*p«»rlc PuU aag^hy PruM. sM ICaunpy, eoMldar. th. 

^reas^ enlaalon. and poaalhl. dagradatln. of aslant air guaUty a- 

rlslng fro. the l^rreeaad ... of coal. Xh. coabnatl.. of coni leada to I 

-ajor amssio*. c£ awlfur and nltrogsn «ddaa. particulate Mtte r. and 

volatile tr«se hazardou. and tooUc asbata.:... Aslant air *i.Uty j 

•taadarda praacrlb. health and welfare snda.«.rl.g level, whan , ,, , * 

the e^t of emaaion. with the attaedit alloueM. calenUted envlzM.' 

■ental lapect for ezpMided coal use will dep«id upon the base Una ewelltv 
«* e*p«.ted gueUty of the air basin involved. 

To estlnat. the air ipiaUty lapect of i«:ree..d coel uee 1, «,t a 
ulnple netter of t3*l«g ratloe of SMk. stack ent«.io.s or air Masure- 
■ents to the .MMt of Increased uee. toosg the f«:tor. thet .wt be 
considered are the type and propertiae of the coal to be weed. Bsetlng 
value. «Ofur. Itrogsn, uih, end trace elsMnt content he coneld- 
ered. The local nl«ro-.t«,rol,gy «wl pin., dlsperalon patterns are In- 
^rtMt in determnlng how Mch of ths poU.ta.ts wUl be found at gronad 
level under nomal as weU as worat-c. co^Utlona. Finally, the loc 

tlon the degree of emaeion. control reglred wUl he deterlaed hr 

both the quality of the area's ..m. 

^ ei' 1^ the potential lnp«-c on 

envlrtsnwiii iilly protected areas In the ^ -'.•'(■atv — 

.11 . 1 - .. / ' -clnlty. Proper preportlonlng 

ell the above factors requires sii Indapth luMwledg. of sMih of the 

wech«>l«» involved as weU as detailed k^wleda. «f the , 

ted whs. coal is bnml. I« e-ltlon. the enrl- bln, steels, 

^t*®e tonic substances forns a sahtle Icme-tMm mmA i _ 

■untie, long-tern and long-raags threat to 

■an and the envlronnant. 

ere. "T ^ T-tlflcatlo. Of 1_ 

cresses la Missions. (2) poesihl. Wt to elr queUty stsnderde. (3) 

l«u*-t«n inp«:ts of tZM. tote et.ta«».. ( 4 ) reseton. of .Of.: dl- 
(S^rgutlc effects to «:ld rain. vlslhUlty degredatlo.. 
chsnge and other onvlronnantt questions. 
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ffi- ffr 10, Air gollimnt Plspwsioa Mod^Uat. by Fahlen, aiMk- 
clzM th* staciw of throo ■•Jor opproKbM lAlch mmf bo used for the 
4 «MUtotlM prodictloo of the of laeteosod cool bomiog oo oir 

i qoollty. Socb colcolocioas rogolro not only o knowledge of theenlselone 

to bo — fron cools of o cortaln chenicol conpcsitlon bnt olso a 
j gnantltotloo netbod to celoto anlsslan catos to sir quality. To do the 

‘ latter ro^iiros the nso of a dlaporalon nodel. an Inportant link between 

MlMlaa ratoa and doao-reaponae or other health effect atodlea. 

Oisporalen nodels ore claaalfled as foUows: (a) Ganssian nodels. 

iililfli ueuMt that the concentration of a pollutant fron a point eowce 
follows the noTMl error cureo; (h) transport nodels. which are based on 
the low of conservation of anas; ar-i (c) stochastic nodels. which are 

based on the laws of probability. 

The Gaussian nodel has very United validity and Its use requires 
knowledge of parsMtera (•‘dispersion coefficients") which cannot be pre- 
accurately. As a result, errors of several bondred percent are 
not iiBirnisi— ■evertheless. It Is widely used and foms the basis for 
the EfA-reconunded •‘aff-the-sbelf" nodels. 

The transport nodels are sore rigorous— especially when cbenlcal 
reactions are Involved— but usually require knowled^ of eddy dlf fusl- 
vlties or "K" values which also cannot be accurately predicted. Since 
I they are leas United than the Gaussian nodels. they are nathenetlcally 

I nore usually require nore conputer tine. Stochastic nodels 

are the noat rigorous and nost adaptable, but these are in a develop- 

T Md require neteorologlcal data that are not always avaU- 

able. The actual selection of a nodel for use in a "preconstruction 
review" Is for practical purpcses llnltad by Ian to gSA-iai inawinilait 
aodsl's. r*— previous ned eli n g studies are discussed. 

conpelllng Is the great need to develop rellAle yet practical 
ntthftdv for the quantitative prediction of the dispersion of the air 
pollutants enltted In coal burning. 


nrmoDOCTioN amd suhmary 
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Chapter U, Atwoaoherlc HodUlcatlona . by Taylor, Chaneldes and 
Green. Is concerned with the atnospherlc Inpact of the release of pollu- 
tants fron coAustlon of fosaU fuels, particularly coal. Such releases 
can cause global perturbations by changing the average conpositlon of 
the earth's lower atnosphere. Hhether such perturbatioos are signifi- 
cant depmds upon the nagnitudes of the conbustlon sources as co^red 
to other natural sources and the rapidity with which these poUutants 
are scavenged before they ate dispersed throughout the atnosphere. 

In this chapter we focus on the potential inpacts of the continued 
or accelerated release of CO^. 110^^. SO^ and aerosols upon global clinate 
and other Inportant envlronnental paraneters. Out present understanding 
of the atnospherlc budgets of these pollutants Indicates that anthropo- 
genic enlsslons of CO^ have already led to «c Increase in global CO^ 
levels, while and S^2 ^ 

While the cllnatlc perturbation inplied by a global CO^ Increase appears 
to be the nost significant global pollution problen we presently face, 
the consequences of Increases in WO^ and SOj upon the envlronnent (i.e.. 
acid rain) are also of concern. In the caiw of atnospherlc particulates, 
the ovIHa tion of S conpounds to produce SO^ aerosols In both the tropo- 
sphere and stratosphere nay ultlnately lead to a significant degradation 
in vlslblUty. Sone alnple techniques for nonitorlng vlslbUlty are 
described. One inportant nethod of nonitorlng aerosols and the deterio- 
ration of visibility on a global scale is the application of renote 
sensing with space technology. This technology is also applicable to 
the global nonitorlng of CO^. SO^. HO^ and the earth's radiation budget. 
The concluding section of the chapter Illustrates this rapidly advancing 
technology. 


26 


ICAAS 


I 


ChapCcc 12, Solid Haata and Trace Elwant li»«ct«. by Bolch, con- 
sldera cbe potential envlrauaental Inp a r ta of tha aolid naatea and txaca 
elenent teleaaea due to coal utilization fot alacttlc ponce* Coal nay 
contain a wide apectcua of tcace elenenta Including Aa, Cd, Ca, Ce, Cu, 
Hg, Ifa. Fb, Sa, Sr, V, Zn, and naturally occnrrlng radioactivity, aape- 
clally the iiraniMn and thoriun aeries. A recent nonograph by Toney re- 
views tha potential lapacts froa trace elenenta, and a conpanlon nono- 
gyapii focuses on the recovery of these waste products aa beneficial re- 
sources. Both nonographa were pid>llshed before the inpact of the le- 
source Conservation and Kecovery Act was reflected in the proposed haz- 
ardous waste regul* Slone. Sone of thf topics not enphaalzed by Torray 
are nore adequately covered in the recent inpact atateaent on the Fuel 
Use Act (DOE, 79) . 

n..r <n g the developaent of the envlronnental regulations of the last 
half of thia decade, there has been an Increasing eapbasls on the leas 
obvious pollutants froa fossil fuel power plants. The return to coal, a 
defined national policy, results la the consideration of the ultiaate 
fate of tcace elements in coal as a significant coal burning Issue. 

This chapter suansrlzes the trace eleaents in coal and fly ash, in- 
cluding the radioactive coaponenta of uranium and thoriun. The fraction- 
ation of coal ash within a typical power plant is presented and the asso- 
ciation of eleaent with various fractions Is discussed. A brief review 
of potential health effects is presented. It is suggested that envlron- 
aental transport and dose-to-rlak aodels he developed and presented In 
oeder to place coal burning on an envlronaental cost scale with ocher 
energy sources. A brief discussion of the Impact of aining, coal clean- 
ing, storage and coal conversion are discussed. Lastly, the importance 
of new laws, naaely. The Toxic Substance Control Act (TOSCA) and Its 
coapanion. The Beaource Recovery and Conservation Act, is discussed. 
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Chapter 13, Aarleulture. by Holtz and Stzeec, describes tbs ssnrlzom- 
ia^act of Increased coal usage on agriculture and natural ecosys- 
ceas. The i^acc of coal residues on agriculcursl and forestry anvlraa- 
aents depends on the partitioning of these aaterlals betweaa bulk solid 
waste (90X) and released emlssian products (lOX). Main cnvironmaatal 
concerns are sulfur ozidss r ule s ss d as atmospheric 90^ sad acid rala fr« 
sulfuric acid aerosols. These subjects have been well documsated quali- 
tatively but not so weU qusntltatlvely. The aagaltude of the offsets Is 
assayed Imdlrsetly throu^ appraisal of visible acute damage to pleats. 
Settlement of specific damage claims for SO^ offsets rsqulrss a eoasldsr- 
atlon of the pertinent features of etiology, environment and degree of 
suaceptlblllcy of cbe plant populations at rlA. 

effects of SO^ are recognlzsd in terns of gaaetic adapta- 
tion cakes place in individual plant spseles in response to airborne 

SOj. Also, the aafceup of populations by species in exposed arses repre- 
sents a biological adapcatliv i At occurs. Changes la genetic and spe- 
cies makeup plant populations any or nay not be deslrabls. Acid rain 
and acid aist from airborne sulfates and alcrates have bean shoM to have 
adverse effects on native and cultivated vegetation, soils, and aquatic 
ecosystems. These effects have increased in recent yssra. Trace elements 
from coal burning nay iapact adversely on agriculture and ecosysteas but 
the effects will most probably be of a much lower order of nsgnicnde than 
Chose of sulfur oxides. 

There Is a potential benefit from disposing of coal solid waste res- 
idue on agcicnltural and forestry soils insofar as the chsmlcal and plya- 
<...1 of Cbe soil are not detrimental to Cbe prodnctlon of 

quality food and fiber. Doe to the high variability of cbe chamical and 
physical nature of coal ash and soil a coapsclblllty anst be ascertained 
before indlscriminatory use of coal ash on soils can be permlttsd. Frs- 
llmlnary studies indicate that certain plant essential elements contaiasd 
in coal ash are readUy available to plants groesi cm ash- e m e nded soils. 
However, there are possible isifavorable c han ges in the soil upon addition 
of coal ash which asst be Investlaged. 
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Chapter 14. Air Pollatlna Health Ettecta . by Schlenker and Jaeger, 
daela with the poaslble health eCfecta of air pollution resulting frow 
large-scale coal utUlzatlan. The pollutants (prlaary ones such as SO 2 . 
particulates. MO 2 . CO. trace eleaents. and secondary ones such as 0^ 
and aerosols) are evaluated according to their physical, chenlcal and 
biological properties. To understand the health Inpacts of these pollu- 
tants pertinent epldenl ologleal studies and controlled laboratory 
studies are presented la which aalaals and hunnaa were exposed to various 
pollutants. Each technique has Inherent advantages and disadvantages. 
Aalnal studies show prlnarlly specles-speclflc responses to exposures, 
hut such studies are Invaluable alrrora of bloaedlcal and aorphologlcal 
changes which nay occur as a result of exposures. For ethical and legal 
reasons, studies using hunan subjects are restricted In length of ex- 
posure and concentration of pollutant. However, the nost Inportaat 
Infomntlon obtained tram such studies Is the acute response by nomal 
and sensitive (swch as astfanatlc) hiuwn beings to well-defined levels 
of single pollutants or conblnatlons of pollutants. Epldeniologlcal 
studies allow one to evaluate the effects of air pollution on large 
medial ■ of people over a llfetlae. Such studies, unfortunately, are 
the nost difficult to conduct. 

I^ldcniologlcal studies can be aidtdlvlded according to the health 
effect conslderad. These sid>dlvlsians are nortality rates, norhldlty 
rates, and cancer rates. Confounding factors such as location and nun- 
ber of nonltorlng sites, anoklng habits of the population (another fom 
of self-induced pollution, elaborated upon In this chapter), geographi- 
cal saasonal variations, socio »%'onanlc factors, occupational exposures, 
and nlgration of Individuals nake Interpretation of epldeniologlcal 
studies difficult. Taking all these factors Into consideration. In nany 
cases the naasnrable relationship of air pollution to health is de- 
creased. There Is evidence, however, that sulfur-dioxide, sulfates and 
particulates, nsjor products of coal coubustion. have sons detrlnental 
effects on the health of children and adults; psrtlcularly sensitive 
subjects react to a larger extejt. 


IMT80DUCTX0N AND SUtBIARY 

Chapter IS. Quantitative Public Policy Asaessnects . by Green and 
Bio, reports the overall approach and principal results of a study 
(ICAAS, 1978) which attenpted to carry out an integrated interdisciplin- 
ary assessnent of air pollution abatenent alternatives In the Tanpa area 
of Florida, a region where coal Is a aejor source of electric power. 

ICAAS *s Initial concept' of public policy decision netfaodologies (PPDH) 
for regulating air pollution began In connection with a proposed Air 
Quality Index (A(^I) Project (KUkAS, 197C). This AQI project was designed 
to be a broad soclo-tecfanlcal research progran leading to the establlah- 
nent of a quantitative scale for air quality. Fron this AQI progran 
plan, we developed our first practical cost/benefit analysis approach 
(ICAAS. 1971). These two overall systens approaches were iaplenented in 
the ICAAS-FSOS s'.-udy (ICAAS-FSOS. 1978) which involved a cha«n of com- 
ponent studlss within the franeuork of two types of public policy deci- 
sion nethodologies on sulfur oxide pollution. One nethodology, the Dis- 
aggregated Benefit/Cost Analysis (OB/GA), was essentially on advanced 
form of economic analysis In iriiich the distributional aspects of S/C 
are considered (l.e. the question of who gets the benefits and who pays 
the costs). The second methodology, the Quantitative Assessment of the 
Level of Risk (QMJt), Is a non-economlc analysis that by-passes nany of 
the difficult problems of translating all important decision factors 
Into monetary terms. 

This chapter concentrates on a few specific but vital facets of 
these system approaches. In particular we discuss quantitative charac- 
terlxatlons of ambient air quality from air pollution health effects 
viewpoint. He also describe a "factor of safety" approach to quantita- 
tive dose response relations based on the use of an air quality Index 
(AQI). These are followed by the application of the QALR-FFDH used in 
the ICAAS-FSOS study which also uses an AQI. He also describe the 
essence of the DB/CA-PFOH used in the ICAAS-FSOS study. Finally, we 
describe some recent dose response results for plants and materials 
whli'h can be used in PFDH's. 
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denand for electric power which Is expected to Increase greatly iron 
1980 to 1995. To aeet this deaand. and also to finance the conversion 
of existing oil and gas fired plants to coal, the utility Industry nust 
raise and Invest unprecedented large suns of aoney. 

At the start of the 1960s, the electric power industry was the 
epltoae of financial strength, with a virtually unlinlted ability to 
raise funds. Today, however, the average coapany la so weak finan- 
cially chat it slaply cannot aeec its capital requlreaents. This de- 
terioration was caused hy a coabinatlon of econoaic and political 
factors — infaltion, both general and in fuel prices, is the root cause 
of the industry's problea. Cost Increases have outstripped productivity 
gains, which has squeezed profits, necessitating rate Increases. How- 
ever, since utility prices are set by regulatory conalsslons, tiae lags 
are inherent iu <d>talnlng rate relief. If costs rise but prices can be 
increased only after a regulatory delay, then obviously profit aargins 
are squeezed, rates of return on invested capital decline, and the coa- 
pany's financial position suffers. 

The situation is aasked by Che fact that Che coapanies now have 
excess capacity chat arose froa Che sudden, sharp reduction in growth 
after 1973. Since this excess capacity has permitted Che coapanies to 
survive and aeet current power demands, Che public has not yet suffered 
to any significant extent. However, excess reserves will soon be used 
up, so that If construction prograas have not been started up well In 
advance, power shortages will follow, accoq>aaled by severs economic 
problems. It is possible, however, to avoid capacity shortfalls. HhaC 
is needed is for utility coaaisslans across Che country Co realisti- 
cally analyze the situation and Chen to allow Che utility coapa ni es to 
charge prices that cover the cost of providing service, including the 
cost of the capital Invested In the plant that provides the service. 


XNTROOUCTIOM AMO SUMMARY 
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Its with 


Chapter 17, Coal and the States; A Puhllc Choice Farsaactiva . 
by Rosenbaua, discusses the Coal Bomlng Issues pciAlea froa the per- 
spective of a political scientist. The devalopaant of a national coal 
policy will depend heavily upon the states for iaplenantatioa. The 
generous discretionary authority exercised by the states in iaplsaantlng 
federal coal policy, together with traditional state powers affecting 
coal use, aesns the states will be aajor actors la any national coal- 
aanagenent prograa. In general, state activities will affect coal util- 
ization through the siting of aiaes, the siting and design of power 
generating plants, the location of coal logistical facilities, the 
enforceaent of air and water pollution staadarda, aad nach else. 

Foiaulatlng state coal policy confronts the state g ov ern — nts with 
competing, and so— tines conflicting, policy objectlv— and policy 
priorities. Aaong the — Jor policy choices that a— t be resolved by 
Che states ace: (a) the relative importance of environ— ntal protection 

anong other policy cAJeetivea in coal n— ; (b) the priority to be 
accorded econoaic growth in coal dev el o p m en t; (c) the distrlbnti— of 
social costs — d benefits fr— coal develop— nt— the "dlstrlbncive 
equities"; and, (d) the rslative priority Co be giv— state or regional 
interests compared to —Clonal on— in cho— Ing c— 1 policy goals. 

Nhile Che federal govern— nt cannot, and should not. attcapt to 
resolve all these — ccers at the state level, federal coal policy c— 
constructively assist the scat— in resolving the— iaauu. bi parti- 
cular, the federal gov e rn — n t should restrain a— sive a— c— 1 utlll— - 
tlon to a few de— d— at a—t, should cr— te a target growth rate figere 
for the 0.S. , should eaph—1— con— rvatl— of e— rgy in the coni — etor 
by dampening demand for n— electric ge— rating f— ilici— and dniuld 
encourage gr— c public involvement in c— 1 policy foranlnti— am— g the 
Western scat— where the environm ent al risks of coal developm— C are 
especially acute. 
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Chapter 18. federal Re ff«latory an«l Legal Aspects, by Little and L. 
Capehart. federal laws which alamltaneously proMte and con- 

strain the Increase In coal use as an alternative fuel to natural gas 
and pctrolew. Congressional desire to encourage or require coal use Is 
evident as early as 1974 with passage of the Energy Supply and Envlron- 

.^1 Coordination Act which set nandatory coal conversion requlrenents 

for power plants with conversion capability. Mote recent laws have ex- 
panded federal authority to require the use of coal by najor fuel 
buralng Ins t a ll atlons* 

Congress has recognlaed that an Increase In coal use will be possi- 
ble only If federal assistance is available to weak links In the coal 
productloo-transportatlon-conboatlon chain. To date, some United 
financial help has been authorised for developing new underground nines, 
for rehabilitating the rail systen. and for buying air pollution control 
equlpMnt. At tUe same tlise. Congress has been standing firm on environ- 
nental laws, trying to protect the health of the pubUc by preserving 
air and water quality standards. Clean air legislation has come In for 
the biggest attack because the control of enlsslons fron coal conbustlon 
places Urge costs on any coal burning facility. The waste disposal 
laws nay also create additional financial burdens If coal wastes are 
teraed hazardous. 

lefon of other fedaral leglsUtlon nay be appropriate if the goal 
of increased coal use is to be net. legulntlon of the transportation 
Industry should be ezaaOned to determine inequities in the Uus which 
favor one mode over conpetlng nodes. Cxanples are unequal federal 
subsidies to various nodes, and varied fomulae for detemlnlng rate 
structures. 

Changes in nlnlag laws have been advocated. The nineral leasing 
program of the federal government has been criticised by almost everyone. 
The mine safety legislation has been accused of decreasing productivity 
. without a proportionate Increase in worker safety. Changes In federal 
law may be necessary to achieve the goal of Increased coal use. but such 
should be carefully examined to determine the side effects as 
well as the expected result. Furthermore, proposed changes should be 
considered In light of InterreUtlng laws and problems. 


CHAPTER 2 


COAL AVAILABIL-'^V AND COAL MINING 

M. Ohaniem and Faremums Fardshiaheh 


I. THE CHARACTERISTICS OF COAL 

Coal deposits have quite variable characteristics, dernding upon 
the original topography and water movenent during Its fornatlon. The 
various types of coal are traditionally classified by a n<a*er of para- 
meters. l nd-«dlng heating value, ash content, moisture content, sulfur 
content, and the division of the organic portion of the coal Into fixed 
carbon and volatile natter. The four main ranks of coal are anthracite, 
bituminous, sub-bituminous, and Ugnlte. Table 1 lists some character- 

istlcs of coal In the United States. 

The American Society for Testing and Materials (ASTM) rank Is based 
on the degree of llthlflcatlon and metamorphism of plant material (DOE 
1979). According to this system, tank is determined primarily by the 
percentage of fixed carbon and the heat value (BTO content) of the coal, 
calculated on a nlneral-matter-free basis. A tm>re detailed discussion 
of the energy quality and other properties of coal based upon the evolu- 
tionary history of coal formation as Interpreted fron an energetic 
analysis perspective is given in Chapter 3. 

As described in Chapter 1. coal la the most abundant non-renewable 
energy resource in the world, and it comprises the major fraction of the 
estimated energy reserves of coal, petroleum and natural gas. The 
United States' share of these recoverable coal reserves is estimated to 
be about 28Z or 4.900 quads (Peters and SchlUing. 1978). (The present 
rate of U.S. energy consumption Is ''•78 quads, of which 18Z Is supplied by 
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